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ABSTRACT 



An analysis was done of central Arctic Ocean acoustic data 
to determine the temporal and spatial character i sti cs of 
transient noise events. Digital ambient noise data from 
the FRAM IV experiment of April 1982 were searched for 
ambient noise transients using a detection program. The 
time series of the resulting detections were examined 
visually to categorize each detection as a transient- 
artifact or false alarm. The transient events were located 
in space using time delays between signal arrival at 
different hydrophones. The cross shape of the FRAM IV 
horizontal array permitted location in both bearing and 
range. The source strength of each event was calculated 
using a simple dipole source model. Refraction and 
scatterring of the acoustic path in the Arctic Ocean was 
taken into account. 



The over a 1 1 number of events detected, and hence their 
interarrival times and spatial density, were all affected 
by the background ambient noise level. The detection 
program used the same threshold signal -tc-noi so level for 
all data tapes, so when ambient noise levels were low more 
detections occurred- The mean interarrival time between 
events was 100 seconds. The interarrival time fit a J 
shaped gamma probability di stri but i on . The number of 
events detected per area decreased with range from the 
array center. Half of the events occurred within 3000 
meters of the array. In this area there were 0.3 events 
per square kilometer per hour. The event population showed 
no predominant angular dependence- The strengths 
calculated using the simple dipole model had a mean of 430 
kM overall and 260 kN during quiet times. Stronger events 
occurred during times with high ambient noise levels. 
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CHAPTER 1 
INTRODUCTION 

Thesis Motivation 

This thesis investigates the spatial distribution, 
strength and rate o-f occurrence at low frequency central 
Arctic noise events- During the last 30 years there has 
been a growing commercial, military and academic interest 
in the Arctic region. 

In this relatively unexplored area there is increasing 
evidence of rich mineral and petroleum resources- Research 
into methods of locating these assets and constructing 
facilities to exploit them have received much attention. 
These facilities must be able to withstand the harsh Arctic 
surroundings. The study of Arctic acoustics helps in 
understanding the Arctic environment and climate- It has 
been shown that there is a direct correlation between 10 to 
20 Hz ambient sound pressure and environmental stresses and 
moments! 101. The ability to use acoustic noise levels as 
an environmental predictor would be a useful tool in the 
protection of commercial Arctic facilities. 

The Arctic ocean serves as a military arena for 
several submarine fleets- The underice environment makes 
detection difficult, increasing the stategic role of those 
fleets- Because of the sound velocity profile of the 
central Arctic there is a surface duct which channels sound 
for long distances- But, the underice profile scatters 
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sound energy, eHecti vely -filtering out high 
•frequenciesCll], The result is that only low -frequency 
signals travel -far in the Arctic, and therefore the low 
-frequency range is the best -for detecting adversary 
submarines- The importance o-f understanding the low 
-frequency ambient noise -field becomes apparent. The actual 
central Arctic ambient noise level is at times much quieter 
than the open ocean, but it contains unpredictable 
transient noise events which inter-fere with conventional 
detection schemes. It has been hypothesized that the 
background ambient noise is the summation o-f these 
transients -from throughout the Arctic basinC73. Analysis 
o-f the spatial and temporal distribution o-f these 
transients is a logical next step in understanding low 
-frequency noise, and improving our submarine detection 
capabi 1 i t y . 

The academic challenge o-f the Arctic lies in the 
sparseness o-f -field data. The Arctic cannot be casually 
sampled. Even simple experiments require expensive 
expeditions. The harsh environment takes its toll on 
researchers and equipment, and reduces the amount o-f usable 
data. Hence, the study o-f the Arctic is like a jigsaw 
puzzle with -few pieces present. The total picture remains 
a stimulating mystery. 

This study analyzes data collected during the FRAM IV 
experiment by Massachusetts Institute o-f Technology and 
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Woods Hole Oceanogr aph i c Institute personnel. The FRAM IV 
ice camp was located in the Barents (Nansen) Abyssal Plain 
at approximately 34° N by 15° E, as shown in Figure 1 — 1- 
The ice was 3 meter thick multi-year pack ice. The ice 
activity was low; there was no ice ridging or lead 
formation around the camp during the experiment. 

The FRAM IV ice camp was set up from 25 March to 11 
May 1982- This study analyzes data taken between March 
27th and April 22nd- The weather was mild, with 
temperature ranged from -35° to —4° C, and wind speed from 
1 to 23 knots. 

The ambient noise was sensed with a large horizontal 
hydrophone array which consisted of two non— uni f orml y 
spaced line arrays, crossing at right angles. The data 
were digitally recorded on a multichannel system- 
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Figure 1-1 Location of the FRAM IV Arctic experiment 
conducted in the spring of 1982. Cl 13 
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Tnesis Contents 

My work in this area began long after the FRAM IV 
ambient noise data were collected. The first step was 
finding the events in the raw data. The data consisted of 
magnetic tapes each containing 20 minutes of digitized 
noise levels. A program was written which searched the 
ambient noise tapes for possible events. Chapter 2 
di scribes the automated and manual techniques used to 
accomplish this detection. 

These events were then located in space using the 
difference in arrival time between hydrophones. This was 
also done with a computer program. The program plotted the 
arrival time delays against range to a trial location, did 
a least squares fit, and chose the location with the best 
fit. This is covered in Chapter 3. 

The peak voltages for each event were used with 
the dipole source model to predict peak source strength. 

The background ambient noise strength was also determined. 
These strength calculations are found in Chapter 4. 

In Chapter 5 the distribution of event interarrival 
time was determined. The event locations and strengths 
were analyzed, and a spatial density found. 

Chapter 6 summarizes the key results of this study. 



Data Collection 
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CHAPTER 2 

DETECTION OF NOISE EVENTS 



Twenty nine FRAM IV ambient noise tapes were searched 
in order to find a population of noise events for this 
study. The specific tapes were chosen from the possible 67 
in order to cover the entire range of days of the FRAM IV 
experiment. However, there were several days when no 
ambient noise tapes were recorded. To help fill these gaps 
five reverberat i on tapes were also searched. These tapes 
were recorded prior to the reverber at i on shot being fired, 
or they were recorded so late in the experiment (SO minutes 
after the shot) that reverber ati ons were no longer present.- 

The FRAM IV experiemnt used the horizontal array of 
omni di rect i onal hydrophones pictured in Figure 2-1. The 
hydrophones were suspended from the ice into the water to a 
depth of 93 meters below the air/ice interface. The two 
crossed lines of the array allowed the possibility of 
localizing events in space. Although 26 hydrophones are 
shown in Figure 2—1, only 24 at a time could be used to 
record data. In most cases a few of the recording channels 
were used for other sensors (geophones or hydrophones used 
in a vertical array) . Most of the time 19 to 21 horizontal 
array hydrophone data were recorded. 

The FRAM IV ambient noise tapes were recorded 
digitally. Figure 2-2 shows a schematic of the system used 
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Figure 2-1 FRAM IV horizontal hydrophone array. C131 
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Figure 2-2 Schematic at the recording system used tor PRAM 
IV data collection. C41 
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tor collecting the data. The input from 
was put through a gain ranging amplifier 
Hz filter. It was sampled at 250 Hz and 
minute magnetic tapeC41. The 24 channel 
dB dynamic rangeC123. 



each hydrophone 
and a low pass 80 
recorded on 20 
recorder had a 120 



Event Detector Program 

The event detector program was written to take digital 

data from a FRAM IV tape and determine where in that tape 
noise events occurred. The program was originally written 
to take data directly from a tape drive, but subsequently 
modified to take the data from a file. The tramread 
program, with a -head switch is used to read the tape into 
the file. This will eliminate any headers and then read 
the digital data straight into a file. A FRAM II tape may 
be read into a similar file using tramread and the switches 
-head and -fram2. The tramread program was written by 
G . Duckworth, and is available to , the Arctic Acoustics 
Program at MIT. 

The event detection programs source codes, flow chart, 
and a short users guide are found in Appendix A. The event 
detection program which reads from a file is called 
hdetect . The detection program reads from a tape drive. 
Both programs are written in the c programming language for 
a UNIX operating system. 

The event detection program follows the block diagram 
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of Figure 2—3. The program initialization portion de-fines 
variables and constants, zeros flags, and requests user 
input such as tape number, date, time and 

channel s (hydrophones) to be used, as well as, the name of 
input and output -files. After this information is 
requested from the user, the program no longer requires 
attent i on . 

The event detection program then reads in a file of 
data, filters the data, squares each data point, and takes 
the square root. The filter was a Par ks-McCl el 1 an digital 
20 to 30 Hz bandpass filter. Its frequency response is 
shown in Figure 2-4. The range of this filter was chosen 
to avoid the Nyquist frequency (125 Hz) and hydrophone 
cable strum (1-20 Hz), and to be compatible with the' analog 
30 Hz low pass filter the data went through before being 
recorded. The data were squared and then rooted to ensure 
positive peak values for ail data points. 

The iext portion of the program used a threshold 
detection scheme to check each channel for possible noise 
events. Far a particular channel a short average of the 
four most recent data points was compared to a long average 
of 64 recent data points. If the ratio of short average to 
long average was over a certain value that channel would be 
flagged for a possible event- The time of the flag and the 
value of the short average were also recorded. All other 
channels averages were taken similarly. If at least 5QX of 



in 
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2—3 Block diagram o-f major madulues a-f the 
ion program- 
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Figure 2-4 Frequency response o-f the Par ks-McCl el 1 an 
digital bandpass -filter used in detection - 
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the channels were -flagged, an event would be declared. 

Then the program would shift to the next time increment of 
data and the process would be repeated. 

A more detailed diagram of the event detection module 
is seen in Figure 2-5. There are four submodules: reset 

flag, set flag, new event, and deactivate old event 
modules- The reset and set flag modules deal with the 
channel flags which trip when a particular channel 
experiences a large si gnal -to— noi se ratio (i-e. the RATIO 
of short average to long average exceeds a certain level) . 
The new event and deactivate old event modules deal with an 
active event matrix which identifies active events, and 
stores channel flag time and amplitude for each declared 
event . 

The reset flag module resets the channel flag if it 
has been more than 0-3 seconds since the channel tripped. 
Spurious peaks on a channel might flag a channel 
prematurely- This reset module prevents a number of 
channels with spurious peaks over a long time period from 
being falsely declared an event. The value of this RESET 
DELAY was determined by examining known events and noting 
that about half the channels tripped within a 0.3 second 
per i od - 

The set flag module determines if the short average to 
long average RATIO has been exceeded and, if it has, the 
module 1) checks to see if the channel flag is already 
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Figure 2—5 Diagram o-f event detection module -Flow and 
decision making. 



trioped, 2) checks to see if this detection is part of a 
recently declared event, and 3) adds its detection 
information to the active event or the channel flag, or 
discards the information depending on the circumstances. 

The short average length of 4 < = 0.016 sec) was chosen so 

that it matched the length of the signal (0.02 sec). This 
provides the maximum signal level since this is long enough 
to get all of the signal and short enough not "p a /erase it 
with the lower surrounding background noise level. The 
long average length of 64 was chosen because a length ratio 
of 15 to 1 had been suggested by KellyC9H for the Large 
Aperture Seizmic Array (LASA) . The detection and 
localization schemes used by this large horizontal array 
were directly applicable to the FRAM IV hydrophone array 
data. The choice of detection RATIO was done through a 
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channel flags is transferred to the active event matrix, 
and the channel flags are cleared. So if a channel has a 
detection and its channel flag is not already tripped, the 
set flag module must first see if the detection belongs to 
a recently declared event. If the active event already has 
that channel flagged, the information is replaced by the 
new detection only if their times differ by less than 0-02 
seconds and the new detection amplitude is greater. This 



means that the existing data can only be replaced by a 
detection of the same signal having a higher peak value. 

If the active event does not have that channel flagged, the 
detection information is entered in the active event 
matrix, and the channel flag is not tripped. The oldest 
active events are checked first, and the detection 
information entered in only one event. If all of the 
active events already have input for tnis channel, anc it 
has been more than 0.02 sec since the most recent input, 
the detection is considered potentially part of an 
undeclared event, and its channel flag is tripped. 

If a detection is made on a channel and its channel 
flag is already tripped, the information in the channel 
flag will be replaced with the new detection only if the 
new detection amplitude is greater. 

The new event module checks to see if at least 50"; of 
the channel flags have been tripped, and if so, declares a 
new event. The 50 X mutual occurrence criterion was used in 
the LASA program with good resultsC91. All the channel 
flag information is transferred to the active event matrix, 
and the channel flags are reset. 

The deactivate old event module was used to remove 
events which were past. A set EVENT DELAY time after an 
event is declared, it is written to the output file and 
erased from the active event matrix. This prevents 



spurious peaks from being added to a event long past. 



The 
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EVENT DELAY time of 0-5 second was chosen because 
inspection of the known events revealed that most channels 
tripped within a 0.3 second period. The RESET DELAY (0.03 
sec) plus the EVENT DELAY result in this 0.3 second look 
for each event. 

The output of the event detection program is an output 
file which contains the time each event was declared, which 
channels were flagged and the time delay and peak amplitude 
for each channel. The time delays are relative to the 
earliest channel flag time, so one of the channels always 
has a zero time delay. 

The first version of the detection program was written 
to take a new short average and long average at every data 
point (every 0.004 seconds) . This program took 4 to 3 
hours to search a 20 minute data tape. A concession to 
speed was made and the program changed to compute averages 
at every fourth sample point (every 0.016 seconds) . This 
reduced the accuracy of the time delays and the ability of 
the program to pick up events- The RATIO had to be lowered 
in order to get the same detections which were obtained 
previ ousl y . 

Studies to find the best si gnal -to-noi se RATIO were 
conducted several times. Development of the LA3A detection 
system had revealed that a 7 dB signal -to-noi se ratio was 
needed for 75/1 detection[3] . This equates to a 5 to 1 
ratio of signal power to noise power. Since I was working 



with pressure vice power I used 2.2 as my starting RATIO. 
This RATIO engulfed me in false alarms. A quick study was 
done around the 2.4 level. The first 10 minutes of tape 
4013 were run at RATIOs of 2.4, 2.45 and 2.48. This tape 
had been visually examined in detail previously, so the 
events were known. The results are shown in Table 2—1. 
Also shown in Table 2—1 are the results of a second study, 
done after the program had been changed to average less 
often. 



Table 2—1 Deter mi nat i an af the Best RATIO 



Average taken 
RATIO 
2.4 
2.45 
2.48 

Average taken 
RATIO 
2.3 
2.38 
2. 4 



data point 
detection rate 
94 7. 

757. 

717 . 

fourth data point 
detection rate 
767. 

717. 

597. 

597. 



false alarm rat 
257. 

187. 

87. 

false alarm rat 
467. 

297. 

337. 

337. 



at every 



at every 



□etection rate = # event detections / # of known events 

FA rate = # non-event detections / total detections 
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Notice how the detection rate has decreased and the 
-false alarm rate increased as a result of only averaging at 
every fourth data point. Averaging less frequently means 
there is a smaller probability that the data paints to be 
averaged will ail lie near the peak amplitude of the 
signal. The signal level is generally lower than that 
detected when averaging every data point, and a 1 ower 
si gnal -to-nci se RATIO must be used to detect the same 
events. But when the si gnal -to-noi se RATIO is lowered the 
false alarm rate increases. 

The RATIO of 2.38 was settled on. This is a 
compromise which gives a detection rate which finds most 
high and medium strength events, and which has a tolerable 
false alarm rate. Because the detection rate is less than 
1007 (717.) there were events present which could be seen 

visually, but were not picked up by the detect iors program. 
The RATIO could have been adjusted to detect all events 
seen visually, but at the cost of a multitude of false 
alarms. The RATIO was kept at 2.38 and used for the 
detection of all data tapes. 

The final version of the hdetect program read digital 
data from a framread file, detected possible events using 
the less frequent averaging scheme, and supplied the event 
time and channel time delays and amplitudes to an output 
file. Once the RATIO had been satisfactorily set the 
program was used to search the FRAN IV ambient noise tapes 
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tor events, and no further program development was done. 
Weaknesses in the program were subsequently discovered, but 
have not been corrected. 

The biggest problem is the accuracy of the event time 
(the time when 50/C of the channels have been flagged). An 
event time reported by the event detection program will not 
exactly match that found by plotting the time series. The 
times are usually within 3 seconds of each other, but have 
been off by as much as 13 seconds in one case. The time 
difference between the two methods is greater at the end of 
a tape, and is likely to occur after a particularly strong 
event has taken place (though there were times when time 
di screpanci es developed without strong events present . and 
also many strong events existed which did not induce 
di screpanci es) . Typically, there might be no time 
difference at the -^irst part of the tape, then after a 
strong event a three second discrepancy would be seen and 
tiiis would be consistent until the end of the tape. 

Because the errors did not appear randomly throughout the 
tape, and because they developed impulsively, I believe 
that the problem lies in the time counter of the event 
detection program becoming offset from the time of the raw 
data, perhaps because of short records in the raw data. 

The event detection program was not written to adjust the 
time keeper in the event of a short record. Correcting 
this may eliminate the time discrepancy problem. 
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It has already been mentioned that the accuracy of the 
time delays deteriorated when the program was changed to 
run more quickly. This also made the detection of the 
event peak amplitude less likely. As a result, in order to 
get accurate locations and source strengths, both time 
delays and peak voltages had to be taken manually from time 
series plots of each event. 

The other major problem of the event detector program 
is that it does not discriminate between an Arctic noise 
transient and an artifact such as an air gun blast or a 
reverberati on shot. Short, strong signals are reported as 
possible events. Adding this di scr i mi nati on to the program 
is the next step in improving its usefulness. 

Vi sual Conf i rmati on 

Visual confirmation was required for all possible 
noise events in order to eliminate artifacts, false alarms 
and multiple detections of the same event. In addition, in 
a few cases visual confirmation revealed two events where 
there had only been one detection. 

The event detection program was designed to preclude 
the need for plotting a time series of each event. The 
output of the program contains time delay and voltage 
amplitude information which can be used directly in the 
location program. However, because of the decreased 
accuracy of the time delay and amplitude information, and 
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because of the event time discrepancy mentioned previously, 
it was necessary to plot the time series of • each event. 

The first step of the visual confirmation is to review 
the tape log for any artifacts that may have occurred 
during the recording. The times are noted, and these are 
compared to the event times given by the detection program. 
Then a time series of the artifact was plotted to determine 
which detections were associated with it. In general, an 
artifact such as an air gun blast did not affect detections 
for over 20 seconds. 

The visual confirmation portion of the procedure 
evolved from a very limited look only at events which could 
not be located with the detection program generated time 
delays, to a three step plotting procedure for each event. 
During the early period of this work the hdetect program 
output was used directly as the input to the location 
program. The location program used the time delays to 
determine the event's location in space. Those events 
which could not be located needed a closer look, and so 
their time series were plotted. The plots were made of a 2 
second period including the event time given by the hdetect 
program. Often there was no apparent event in this time 
series plot, and the detection was declared a false alarm. 
When an event did plot, manual time delays were taken and 
used to locate the event. These manual delays located 
these events with better accuracy then the hdetect 



generated time delays. 



It soon became apparent that the 



best answers would be obtained by taking manual time delays 
of al 1 events. Trying to localize the events with the 
program generated time delays was dropped from the 
procedure, and the first step after getting the hdetect 
program output became doing a 2 second time series plot of 
each event. 

After a dozen tapes had been analyzed in this manner 

the discovery of the event time discrepancy was made. 

After plotting a dozen events right at the time shown by 
the program, the final two dozen event of tape 4009 all 
appeared to be false alarms. The quality of the tape was 
good (low background noise), so this seemed highly 
suspicious. A broader search of the time around each event 
showed that the final two dozen events were not false 
alarms, but were events with times 5 seconds different than 
those indicated by the program, so that none of those 
events had shown up on the 2 second time series plots. The 
method of visually confirming events was changed so that a 
waterfall plot was made around the time of each event. 

This showed the exact time of the event, and helped discern 
the pattern of time discrepancy between time series and the 
hdetect program. Once the pattern was found events were 
easy to find and false alarms could be noted. The 
confirmed events were then plotted with 2 second time 
ser i es . 



The -final step o-f visual confirmation was simply 
plotting and replotting the events to the proper gain so 
that the higher voltage amplitudes would not be cut off. 
These peak voltages were manually taken -from the time 
series plots and used to determine source strength. 

The -final method used -for visual conf irmation was: 

1) Check tape log -for artifacts, and eliminate 
those from further analysis. 

2) Plot a waterfall time series around each 
passible event, separate real events from false alarms, and 
find true event times. 

3) Plot 2 second time series of each real event, 
adjusting gain to keep from clipping higher voltages. 

Using this technique certainly reduced the false alarm 
rate. A breakdown of the detection statistics for tapes 
that had been examined by both methods is found in 
Table 2-2. 

Table 2-2 Detection Statistics for Two Visual 
Conf i rmat i on Methods 

Artifacts Events False Alarms 

□ rignai Method 16.97. 37.57 45.67 

w/ 2 sec plots 

Ultimate Method 16.27 65.47 18.47 

w / w at er f all p 1 o t 
2 sec c 1 ot 

(Percentages of detections classed in each category) 



The human interpreter was a necessary tool in this 
scheme- There was not necessarily a one-to-one 
correspondence between events and detections. There were 
cases where a strong event would cause multiple detections, 
and cases where two events occurred at the same time and 
caused only a single detection. In some cases a series of 
detections seemed to be an event and an echo, or perhaps 
straf . This would be counted as a single event. 

The method of determining whether a detection was a 
false alarm or a weak event was sometimes difficult- In 
general, if the detection program indicated a possible 
event, "something" could be seen on the waterfall plot. 

The detection was dismissed as a false alarm if no pattern 
for taking time delays could be seen. (Because of the 
shape of the hydrophone array there were consistent 
patterns of time delays depending upon the direction to the 
event.) Presumably the false alarm rate depends upon the 
training and attention of the human interpreter. 

Manual time delays were taken from an arbitrary 
reference to the crossing of the largest peak to peak 
amplitudes, as shown in Figure 2-6. For most events this 
was clear, but for weak or complex events some intuition 
was needed. 

Voltage amplitudes were taken as the maximum peak 
voltage in the event signal. All were taken as magnitudes 



regardless of sign. 



-time delay 

e —i 




Figure 2-6 Sketch showing point of measurement for event 
time delays. The event is timed at its zero crossing 
between the largest pair of positive and negative peaks. 
(The measurement/analysis system has a polarity of negative 
voltage for positive pressure. ) 



Ail types of noise event signatures previously 
observed by DverC7U were seen in the ambient noise tapes I 
evaluated. The majority of events were pops and extended 
pops. There were also a few whines and straf events. 

While signature types were noted in general, the signature 
type of each individual event was not recorded- 



CHAPTER 3 



LOCATION OF NOISE EVENTS 

Event Location Program 

The program used -for localization was based on the 
program FQUAK by Peter Stein C131- This program places the 
event at different trial locations and computes the slant 
range to each hydrophone. Figure 3-1 shows the coordinate 
system used for these cal cul at i ons . These are plotted 

against the time delays and a least square fit is done to 
determine slope as shown in Figure 3-2. 
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The standard deviation of the time delays from the slope 
line is figured. 



si gma 




CT 



(3-3) 



The location having the lowest standard deviation is the 
location of the event. The inverse of the slope is the 
group speed of the signal. The y intercept of the plot is 

added to the reference time of the manual time delays to 
get the time the event actually occurred (as opposed to 
when it reached the hydrophone array). 
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Figure 3-1 Coordinate system used for calculation o-f event 
1 ocati on . 



Time Delay (seconds) 
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Figure 3—2 Least, squares -fit o-f time delays and test 
location slant range. 
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I have assumed that the strongest peak pressure sensed 
at the hydrophone is due to a waterborne acoustic 
propagation path -from a source located in the ice sheet. 

The signal enters the water near the source and propagates 
directly toward the hydrophone. 

I have assumed that the signal does not bounce o-f-f the 
ocean bottom or the ice canopy before reaching the 
hydrophone. Paths bouncing o-f-f the bottom would produce 
signals with much 1 ower energy than the direct path signal, 
and can be ignored. Signals bouncing off the ice canopy 
are too energetic to ignore but, as I show subsequently, 
they do not affect the time delay computations significantly. 

The location program is based on arrival times being 
related to slant range, R , and does not take the upward 
refraction of the acoustic path into account. The impact 
this has on the results is discussed in the next section. 

The location program takes as input a file of time 
delays and voltage amplitudes, and outputs a file 
containing the best event location, sound speed, and 
standard deviation. It also computes source strength based 
on the voltage amplitude inputs, the event location and a 
spherical spreading loss. This feature was originally 
included so that the source strength could be computed 
directly from the event detection program outputs- Since 
the peak voltages recorded by the detection program are not 
as accuraLe as those done by hand, and since the 
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transmission loss does not -follow simple spherical 
spreading, these computed source strengths were not used 
-for any part o-f this study. 

The FQUAK program set up a grid of points around a 
specified center position. The grid consisted of a point 
every 100 meters from —5000 to +5000 meters in both the x 
and y directions. This resulted in 100 x 100 test 
locations. When the best test location was found the 
interval spacing was reduced to every 10 meters, and 
another 10,000 test locations were generated using the best 
location of the first round as the new center. The process 
was repeated with a 1 meter interval to get the final 
answer % The scheme evaluated a total of 30,000 test 
locations, covered a range out to 5000 meters, and took 
about 20 minutes to run. 

I noted that a significant number of events found with 
FQUAK were at the range limit of 5000 meters. The program 
location was written to search a larger area faster. The 
fineness of the grid was decreased to 20 x 20 vice 100 x 
100. A 1000 meter interval was added to enable the program 
to search out to 10,000 meters. This reduced the total 
number of test locations to 1600 (20 x 20 x4) , and the time 

to one minute. location gave answers which were very 
consistent with FQUAK - except in one particular situation. 

The wider grid size led to one problem. The location 
program sometimes found the lowest standard deviation for a 
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point in the quadrant directly opposite the true location. 
This was suggested by the sound speed being reported as 
approximately -1440 m/sec, as illustrated in Figure 3—3. 
This problem was solved by modifying location to make the 
program f ine locate . This program used the grid size and 
spacing of FQUAK, and centered the search so that the user 
could designate which of the quadrants would be searched. 

A casual leak at the manual time delays of an event easily 
reveals the appropriate quadrant. This program works well, 
but is as slow as the original FQUAK* It was used rarely. 

As with the original FQUAK , I began to notice that 
some events were located at the range limit of the location 
program. This led to the modification of the location 
program to form the program far locate . This program uses 
the location grid size and fineness, but allows the center 
point to be any of the far corners of the original location 
grid, or at the limit range at each of the cardinal points. 
This is shown in Figure 3-4. This allowed events to be 
located out to 20,000 meters - 

The location program source code and a brief user's 
manual are found in Appendix B. This program was written 
in the c programming language for the UNIX operating 
system. This program was developed to the point of 
usefulness, and then used to locate events. No further 
program development was done (except the very minor changes 
to produce finelocate and far locate ), so there are surely 
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3—3 Least squares -fit to a false location 180° 



Fi gure 
away. 
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locat ior: 



f in e locate 




far locate 

Figure 3—4 Areas covered by location - far locate and 
f ine locate . 
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improvements to be made. 

The location program is quite interactive. One 
hydrophone with a bad time delay can change the slope and 
location a great deal. An event is located by eliminating 
bad time delays and checking the sound speed and standard 
deviation of the location. In some cases no hydrophones 
needed to be removed, but in most cases at least one 
hydrophone was removed before an event was considered 
located. The sound speed was the major indicator of 
whether an event had been located. If the sound speed was 
between 1380 and 1500 m/s the event was considered located, 
□f course an attempt was made to get close to 1440 m/s ec . 
This had to be balanced with reducing- the standard 
deviation. A standard deviation below 0.01 seconds was 
considered good. 

A table summarizing all of the events and their 
location parameters is found in Appendix C. The standard 
deviations (sigma) are given in two sets of units. The 
first is the sigma calculated by the location program, and 
it is in seconds. The second sigma is a translation of 
that standard deviation to meters using the sound speed 
calculated for each particular event. The standard 
deviations ranged from 0.0010 to 0.0327 sec, with 0.0077 
sec being the average- The significance of this standard 
deviation will be discussed in the next section. 



In some cases just removing suspect time delays did 
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not lead to a localization- A reexamination of the event 
time series was done to see i f any of the manual time 
delays was incorrect- Often a reexamination of the time 
series produced a change of 1 to 4 of the time delays. 
These corrected values plus values from the other channels 
would then be used to locate the event- About 15X of the 
events required reexamination. Most of those were 
subsequently 1 ocated - 

Despite the above efforts, there were a few events 
that could not be located within the 1380 to 1500 m/sec 
sound speed limits- These events may be from propagation 
paths other than the assumed direct acoustic path. Events 
arriving primarily through the ice longitudinal wave or 
the ice flexural wave would have phase speeds above and 
below my sound speed limits- These non— 1 ocat abl e evvents 
are indicated in the event location summary of Appendix C, 
and they were not used for any analysis which required 
accurate 1 ocat i on - 

Figure 3—5 shows the position of the events located 
within a 2 km square centered on the array origin- Figure 
3-6 shows the position of all events located- 
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Figure 3-5 Noise events located within a 2 km 
surrounding the array origin- 
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Figure 3 6 Position o-f all noise events located. 
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Effects of Refraction on Location 

A sound speed profile was used to get the refractive 
paths for various ranges- This was simplified by the fact 
that all the hydrophones were at a depth of 93 meters. 
Assuming that only the “direct” path is involved means that 
each horizontal range has to have a unique launch angle in 
order to reach the hydrophone at its specific depth. Rays 
were launched into the layers of the sound velocity profile 
and the horizontal range to the hydrophone was calculated. 
The time required to travel the refractive path can be 
calculated and compared to that of the slant range- This 
time error can then be related to the error of the location 
program. 

Figure 3-7 shows the linearized sound velocity profile 
that was used. It is based on the sound velocity profile 
reported for the eastern Arctic ocean by ChenCll. Figure 
3-3 helps to illustrate the scheme used to calculate the 
ray paths. Equations 3.4, 3-5 and 3.6 were used to 
calculate angles, ranges, depths, and propagation time. 
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1.434 1.438 1.442 1.446 1.45 1.454 1 . 4-58 1.462 1.466 

(Thousands) 

Sound Speed (meters/sec) 



Figure 3-7 Sound velocity profile used in predicting 
refractive paths. 
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Figure 3—8 Scheme used to compute retractive propagation 
paths. Source is located at z = Om and hydrophone at 
93m. Sound speed gradients change at 80m, 254m and 36 
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Llsing the known depths (z) and estimated speed 
gradients (g) of Figure 3-7, and choosing a particular 
launch angle <S Q ) , all of the subsequent angles of 
intersection of the layer interfaces ( 0 ^ , Q-- , 8 ^) can be 
found from Equation 3-4. The angle which intercepts the 
hydrophone at its depth ( 82 ) can also be found- Uiith the 
angles known, the horizontal range and propagation time for 
each layer can be determined with Equations 3—5 and 3—6. 
These are combined to get the total horizontal range and 
propagation time. It should be noted that there are two 
ranges and times for each launch angle. The first path is 
that which intercepts the hydrophone on the way down, while 
the other intercepts the hydrophone as it is refracted back 
toward the surface. The maximum depth reached by the 
propagation path was also found, and those paths that went 
below 754 meters, resulting in a range greater than 30600 
meters were not reported. A tabular summary was made of 
launch angle, horizontal range from the hydrophone, maximum 
depth and propagation time, and this may be found in 
Appendix D, along with more detailed tables listing 0^_^, 

r 1-5 and t l-4' 

Rays connecting source and hydrophone with one or more 
bounces from the ice were not considered here. The effect 
of those rays will be taken into account in Chapter 4. 

The refractive propagation time was less than the 
slant range propagation time because the refracted path 
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travels through faster water- The slant range propagation 
time was calculated by dividing the slant range by 1438- 48 
m/sec - the average sound speed between 0 and 93 meters 
depth. The time difference between the slant range path 
and the refracted path are shown in Figure 3-9 as a 
function of horizontal range. 

This time difference is greater than the average 
standard deviation of the location program only after 
13,000 meters, and the time difference at 20,000 meters is 
only about twice that average. The standard deviation does 
not reflect the time difference due to refraction because 
all of the hydrophone time, delays are adjusted in the same 
manner and direction. Figure 3—10 shows that sigma does 
not grow with horizontal range. Refraction effects do not 
influence the standard deviation greatly. Closer than 
13,000 meters the range errc-r caused by other factors masks 
any error from ignoring refraction. 

There is a better point of focus for examining the 
effect of refraction, and that is the change in time delay, 
not the change in the propagation time itself. A point was 
chosen at approximately 5 km from the origin of the 
hydrophone array, and another chosen at approx i matel y 6 km. 
The slant range propagation times and refractive 
propagation times were calculated for each point. The time 
delay between these two points was 0.6857 sec for the 
refractive path and 0.6930 sec for the slant range path. 
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Figure 3-9 Time difference between the slant range 
acoustic path and the refractive path as a function of 
horizontal range. 
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Figure 3—10 Scatter plot of standard deviation and range 
ior all located noise events. 
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The difference was 0-0073 (about IX), or a deviation of 
0.0037 sec -for each hydrophone. Two points at 
approx i matel y 18 and 20 km were also evaluated. The 
di-f-ference between their time delays was 0.0182 sec (about 
0. IX) , or 0.0091 sec per hydrophone. These numbers are the 
same order as the total error o-f the location program. 

The location program may compensate -for some o-f this 
error by raising the sound speed. If just the points above 
were used, the sound speed would go from 1438.5 to 1453.8 
at 5 km and to 1458.6 at 20 km. With 24 time delays being 
used in the location program the effect may not be as 
great . 

The main source of error in the location program is 
the quality of the manual time delays. When the signal— to- 
noise ratio was low, picking the correct peak was often 
difficult. The standard deviation will reflect the 
judgement of the person picking off the time delays. The 
time delays were only measured to the closest 0.003 sec. 

It is interesting to note that 0.006 seconds equates to the 
width of a pencil tip on the time series plot scale. 

The final question to be answered is "How do the 
standard deviation and refraction errors equate to the 
range and bearing accuracy of the location program. " Two 
hypothetical noise events were investigated, one at 5000 
meters (2845, 4136) and the other at 20,000 meters 
(-8253, 18896) The time delays for slant range propagation 



and refractive propagation were calculated- The location 
program was run for each set of time delays, and for each 
set partially contaminated with 0-016 sec errors- (Zero, 
+0-016 and —0.016 were each added to one— third of the time 
delays.) The results are summarized in Table 3—1- 

Table 3-1 Results of location Program Accuracy Test 
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The refraction contaminated by errors case is closest 
to what was input into the location program for the field 
events. This table gives an estimate of the accuracy of 
the location program as 800 m at 5 km, and 2000 m at 20 km. 



The bearing accuracy is excellent- 
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CHAPTER 4 



STRENGTH OF NOISE EVENTS 



Acoustic Source Model 

The dipole is considered a possible source model. 

Peak values -for the source parameter of -force, F , are 
used . 

The acoustic pressure due to a non— convect i ng compact 
dipole source, in a nonr ef r act i ng infinite medium, isL5I: 



Figure 4-1 shows the orientation of the presumed dipole. 

The angle B is the launch angle from the horizontal plane 
down into the water. 

Assuming that F may be expressed as a harmonic, 




(4-1) 



where R = slant range. 




the pressure may then be expressed 



as 



sin 0 27T-FF F 



(4-2) 
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Solving -for F gives 
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In the -far field the 2'TTfR term dominates the sum in 
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Figure 4—1 
de-f i ni tion 



' / 




Assumed dipole source orientation and 
o-f launch angle, Q - 
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the denominator- For the lowest -Frequency considered in 
this study (20 Hz), c (1440 m/s ec) is 10/1 at 2 TT t F; at 115 
meters and only 1-17- o-f 27TtR at 1000 meters. When the c 
in the denominator is neglected the -force can be written 
as; 



2P 0 R\ 

sin 8 



(4-4) 

where is the wavelength. 



The peak pressure, p Q , should lead to the peak force, 
F q . This definition of force was used as the parameter 
for dipole strength- Event signatures that were recorded 
from a source within 300 m of a hydrophone were not used to 
calculate dipole strength, F Q , from peak pressure, p Q - 
For this model the peak acoustic pressure must be 
iound. The hydrophone sensitivity of —159 dB re 1 volt per 
1 fJL Pa was used to convert voltage to pressureC 171. 



1 volt => 89 N/m £ = 89 Pa 



i n — S'. 



The dipole strength formula requires wavelength, \ 
Frequency was taken from the time series plots for each 
event via axis crossing rate, and \ was determined by 
dividing c (1440 m/sec) by the frequency. 



Launch angle needed for the dipole model can be founn 
as in Chapter 3 by assuming a sound velocity profile and 
computing the refractive path. There is a unique launch 
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angle for each horizontal range when the path is purely 
refractive, but a range of launch angles when surface 
reflection paths are included. 

The final parameter in the dipole strength formula is 
slant range. The answers obtained using slant range are 
the strengths based on spherical spreading in a non — 
absorptive medium, equation 4-4. Because the spherical 
spreading assumption is a poor one, (refractive and surface 
reflective propagation paths are caused by the Arctic sound 
velocity profile), equation 4—4 must be modified. The 
effect of refraction on spreading loss will be discussed in 
the next section. 

Volumetric absorption was found by using the 
absorption formulas of DyerCSD. Assuming a pH of 3.2, a 
salinity of 33.5 °/oo, a temperature of 0° C and a pressure 
of 40 atmospheres , I calculated the total volumetric 
absorption to be 1.3 x 1 0 dB/km for an 80 Hz signal. For 
my maximum horizontal range of 20 km, the absorption would 
be 0.026 dB. This is not significant, and I therefore did 
not include a volumetric absorption correction in the 
strength cal cul ati ons. 

Effects of Refraction on Transmission Loss 

Spherical spreading loss in a nonref r act i ng medium is 
illustrated in Figure 4— 2C14I1. The sound pressure is 
pressumed to spread radially. The sound pressure squared 
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Figure 4-2 Spherical spreading. 



(From Ur ickC143) 
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is proportional to intensity, and intensity is the power 
per unit area. Since the power -from a source is constant 

P = I 1 47TR 1 2 = I 2 47TR 2 2 - (4-6) 

The intensity at the reference range of 1 meter can be 
related to other intensities by 

*ref 

I c . = ~ . (4—7) 

47 T R* R^ 

O 

Since I = p^/y0c , this can be expressed in terms of 
transmission loss, H . 

P R^ - 

H = -10 log = 10 log =20 log R , (4-6) 

D r ^ * 

r eT 

in dB r_e the distance reference, taken as 1 m. 

The spreading scheme for a refractive medium is shown 
in Figure 4— 3C23. This is based on ray theory which 
assumes that acoustic energy does not cross the rays, with 
energy contained between two rays being conserved. The 
intensity at the reference range between the two rays shown 
is: 



I 



P 



P 



27 TRCQ 5 0 O R^S 277* cos a Q ^/\a 



(4-9) 



At a horizontal distance r meters from the source, the 



intensity is: 
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1 ) 



Spreading loss in refraction. 



("From Clay and 



T = 
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(4-10) 



The relation between intensities becomes: 
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The loss due to spreading is: 
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In terms ot pressure, I r = p r ^/p 1 c 1 , and I ref = P re +^/ p o c a 
and therefore 
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since Pi'Po = 1 in seawater to an 

excel 1 ent approx i mat i on . 



Applying the dipole model to this spreading loss 
equation gives 



in “(3 
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Equation 4—15 assumes a unique refractive path between 
source and hydrophone. In the Arctic there may be other 
paths due to the non-specular scattering o-f rays off the 
ice canopy. Figure 4-4 illustrates how rays normally 
trapped in a surface duct, may be deflected down to a 
hydrophone. The minimum vertexing angle calculated -From 
the linearized sound velocity profile of Chapter 3 was 
0.064 radians. For a ray to stay in a surface duct above 
the hydrophone it must be reflected from a slope of less 
than 0.032 radians or about 2°. It is reasonable to assume 
that the ice canopy lacks local levelness to this order, so 
that non— specular rays must be accounted for. 

The rays which rebound from the ice canopy experience 
some loss. The attenuation for the FR’AM IV experiment has 
been reported at 0.1 dB/km at 80 Hz E 1 1 1 . This attentuation 
may be converted to a loss per bounce. 



The cycle distance depends on the launch angle and the 
sound speed gradient. For a launch angle of 0.032 radians 
and the assumed sound velocity profile of Chapter 3, the 
cycle distance is 3.7 km. Therefore, the loss per bounce 
is about 0.4 dB. This loss is low enough that even a ray 




b 



dB/km 



( 4 - 16 ) 



X 



where b = loss per bounce, and 
X = cycle distance. 
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Figure 4-4 Specular and nan-specular reflections from the 
i ce/ water interface. 






which has bounced several times may contribute a 
significant amount of energy at the hydrophone. The non- 
specular rays cannot be ignored. Bounce loss at 
frequencies less than 80 Hz are even smaller, since the 
data show a roughly linear dependence on frequency. 

To account for the non-specular rays the spreading 
loss is calculated using the ray averaging techni queCSl . 
The pressure from a particular ray at a given depth and 
horizontal range, assuming a dipole source model, is 



sin^Q 



p^(r,z) = 



o d8 dr 



r I tan Sj | dr X/2 



(4-17) 



dr 

The term represents the probability that a rav bundle 

X/2 

will cross a certain depth, as shown in Figure 4—5. For a 
single linear sound speed gradient the cycle distance can 
be written as 



2c 

X = 1 sin Q = 2r sin Q_ * 2^3 , (4-18) 

g O C O CO’ 

where r c is the radius of 
curvature, to a good approximation constant 
for all small angle rays in a linear sound 
speed gradient. 

Applying equation 4—18 to equation 4—17, and using the 
small angle approximation, gives 



p^(r,z) = 



S„ d8. 



r_ r 6 . 



(4-19) 
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Figure 4-5 Probability of a ray bundle crossing a certain 
depth at a given horizontal range. 
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In order to average the contributions o-f the possible 
rays- this pressure is integrated over all possible angles 
for a given receiver depth, and then averaged over depth 
down to the hydrophone at z Q . 




where 0 V is the maximum launch 
angle o-F a ray that will hit the hydrophone 
at a given range, and 0 m is the minimum 
launch angle. 

The angle 0^ is a -function o-f 0 Q and z 

z = r c Ceos 0^ - cos 0 Q H . (4-21) 

Using the small angle appr ox i mat i on -for cosine leads to 




Substituting this into equation 4-20 and evaluating the 
integral over angle gives 




= — for all z so the second term within the integral 
is always zero. Evaluating the first term over depth gives 
an expression for pressure in terms of r and 8 V - 
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(r ) = 



2 A^ 



-> r z. 



0 , 



[ 2z 1 ~ l/ j:_ 

- d 



(4-24) 



This expression can be used to -find the source strenath. 



r? r> 

A- = p- 



o r z , 





o, -> 3/2" 




—^ol 


8w" “ 


q r - — 


V 


v r_ 




c J 



= G (r ) p ^ 



(4-25) 



The spreading function, G, is presented as a function of r 
alone since Q v depends on r. For each r there is a unique 
Q v , and therefore, a unique G. The spreading function was 
calculated for horizontal ranges from 300 m to 20,000 m, 
and tablulated in Appendix D. ‘The spreading function is 
shown in a log— log plot in Figure 4—6- For comparison the 
equivalent spherical spreading for a dipole source is also 
shown. From this one can see that source strengths 
calculated using the spherical spreading law lead to an 
unrealistic dependence on range. 

In order to get 8 V/ and r Q a linear sound speed 
gradient of 0-054 sec ^ was chosen- This gradient gives 
the same 0 V/ at r = 3 km as the multiple step profile used 
in Chapter 3. Three kilometers was chosen since it was the 
median horizontal range for the noise events. 

The spreading loss function and measured peak pressure 
magnitudes were used to calculate diode strenath. 



F q = 2\A = 2 \ p 0 ^G (r ) . 



(4-26) 
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Figure 4-6 
1 oss - G (r > , 



Spherical (heavy line) and refractive spreading 
as a -function of horizontal range. 
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For a particular event the source strength was 
calculated from each hydrophone, and then the average taken 
as the source strength for the event. The standard 
deviation within each event ran from 97 to 1257 of the mean 
value. An event strength summary in Appendix C lists the 
mean measured peak pressure and the mean dipole strength, 
along with the standard deviation on these values for each 
event . 

Strength of Background Noise 

I was interested in the effect that envi r onmental 
loading might have had on the temporal , spatial and 
strength statistics. It has been shown by liakris and 
Dyer C 101 that low frequency <10—20 Hz band) ambient noise 
rms pressure, averaged over a long time, correlates well 
with environmental stresses and moments. Since I had 
ambient noise pressure for most of the period of the FRAM 
IV experiment, and since I had environmental stresses and 
moments available for only a part of the time, I chose to 
use the 20—80 Hz 1 ong— t i me— aver age rms pressure as my 
envi r onmental i ndi cat or . 

The 10 to 20 Hz band ambient noise pressure was 
converted to 20 to 80 Hz band pressure in the following 
manner. Figure 4-7 shows the typical spectrum for central 
Arctic pack ice noise. The portion of the spectrum between 
10 and 100 Hz can be approximated by a straight line. 



SPECTRAL DENSITY, Po*Hr‘ 
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Fig are 4-7 Composite central Arctic ambient noise spectrum 
observed during the FRAM IV experiment. (From DyerC61) 



log S = A C log f ] + B , (4-27) 

where A = slope = -1.7273 Pa^/Hz^ , 

B = intercept = —1.0909 Pa^/Hz , 



or 



S = 10 B f A 



(4-23) 



The band rms nressure relates to the soectral level by: 



P rms,b 



/ 



S df 



(4-29) 



I have assumed that as the sound pressure level changes 
-from time to time the intercept fc B changes, but the slope 
remains the same. By substituting equation 4-28 into 4—29, 
and using the known 10 to 20 Hz ambient noise band, B can 
be written in terms of the known oressure. 



B = 1 QG 



P^rms . 10-20 



(4-30) 



where 



K i = 



20 (A + 1> - io (A + 



A + 1 



= 0. 1020 Hz 



The ambient noise r ms pressure for the 20 to 80 Hz band may 
now be found. 



80 



P rms. 20-80 



- / 



10 B f A d-f 



(4-31 ) 
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k 2 

rms, 10-20 T7~ ? 

K 1 



where 



go (A + 1) 20 (A + 1) 

A + 1 



0.0988 Hz , 



or finally. 



p rms, 20-80 



P 




(4-32) 



p rms, 20-80 ~ 0.98 P rmS5 io-20 

Thus the band from 20 to 80 Hz is virtually identical to 
the one from 10 to 20 Hz in rms pressure, for long-time- 
averages, and in turn, is an acceptable surrogate for 
environmental forcing (applied stresses and moments). The 
20 to SO Hz band ambient noise rms pressure for each of the 
tapes investigated is found in Table 4—1. 
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Table 4-1 20 to 80 Hz Band Ambient Noise rms Pressure 



Tape # Date Recorded Prms ^0-80 

’(Pa) 



4001 


3-27-82 


Not Available 


2001 


3-29-82 


Not Available 


2009 


3-30-82 


0. 022 


3001 


3-31-82 


0.019 


4003 


4-01-82 


0 . 044 


4005 


N 

CD 

1 

rH 

0 

1 


0. 035 


4007 


4-01-82 


0 . 022 


4009 


4-02-82 


0 .010 


4011 


4-02-82 


0.010 


4013 


4-03-82 


0.013 


2023 


4-08-82 


0.037 


4015 


4-09-82 


0 . 040 


3047 


4-13-82 


0.010 


4016 


4-15-82 


0.017 


4019 


4-15-82 


' 0 .016 


4021 


4-19-82 


0 .013 


4023 


4-19-82 


0. 01 1 


4024 


4-19-82 


0 .011 


4027 


4-20-82 


0.012 


4029 


4-20-82 


0 .012 


4031 


4-20-82 


0 . 0 1 2 


4033 


4-20-82 


0.012 


4040 


4-21-82 


0 . 034 


4047 


4-21-82 


0. 1 14 


4049 


4-21-82 


0. 140 


4051 


4-21-82 


0. 140 


4053 


4-22-82 


o . oao 


4055 


4_22— 


0.032 


4057 


4-22-82 


0. 053 


4059 


4—22—82 


0 . 065 


4061 


4-22-82 


0. 034 


4063 


^ Q'? 


0.023 


4065 


'■?'“> 0 ,— k 


0.027 


4067 


^ 0 *7? 


0 .013 
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CHAPTER 5 

ANALYSIS OF NOISE EVENTS 
Detection Anal ysi s 

A total of 34 tapes was examined, for a total time of 
662 minutes- (For a few of these tapes the entire 20 
minutes was not used.) 

There was a total of 499 detections of events flagged 
on at least 50"/- of the hydrophone channels- Of these, 139 
were man-made artifacts, and 125 were false alarms 
(detections which were so weak that no pattern for taking 
time delays could be discerned). There were 199 unique 
events, and 36 multiple dtections of those events. Stated 
in another way, of the detections which were not artifacts, 
65- 3X were strong enough to support analysis and 34-7"/ were 
too weak to reasonably analyze, and hence labeled false 
al arms - 

Since the detection process depends on si gnal — t o-noi se 
ratio, the level of background ambient noise should affect 
the event detection rate- Figure 5-1 shows normalized 
ambient noise pressure, number of false alarms per tape, 
and number of unique events per tape for each tape 
examined- There is some trend for more events being found 
when the ambient pressure is Ion, and more false alarms 
declared Nhen the ambient pressure is high * 

This is more clearly seen in Figure 5—2, which shows 
the average number of false alarms and unique events found 
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Julian dot* 

4 

Figure 5-1 Normalized ambient noise pressure, number of 
unique events per tape, and number of false alarms per 
tape for each data tape examined. 
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1771 Events 



rms Pressure (P Q ) 

False Alarms 



Figure 5-2 Average number o-f -false alarms and uniaue 
events per data tape -for -four ranges o-f ambient noise 
pressure . 
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per tape in each o-f four background noise pressure ranges. 
The 0.01-0.02 Pa range used 15 tapes to compute its 
average, the 0.02-0.03 Pa range 4 tapes, the 0.03-0.04 Pa 
range 5 tapes, and the over 0.04 Pa range 8 tapes. Two 
tapes were recorded during the first few days of the FRAM 
IV experiemnt, before the 10-20 Hz band ambient noise 
recordings were started. 

A breakdown of detections for each tape is found in 
Appendix C. 

Temporal Analysis 

The inter arrival time between events ranged from 1 to 
1064 seconds. Each event time was taken to the nearest 
second, and no events were taken as having the same event 
time. If two events happened in the same second, one was 
judged to be earlier, and the two events were given event 
times one second apart. The interarrival time for a 
particular event was measured from the previous event, 
except for the first event of a tape, which was measured 
from the start of the tape. 

The interarriva! times were divided into bins of 20 
seconds. The first bin ("0") contained events which had 
interarrival times from 0 to 19 seconds, the second bin 
from 20 to 39 seconds, and so on. The number of events per 
bin is presented in Table 5-1 and shown graphically in 
Figure 5-3. A complete listing of interarrival times for 



Number of Events 
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Figure 5—3 Number o-f events -found per interarrival time 

bin. 
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each event is -found in Appendix C. The mean of the 
interarrival times ( JJi ) is 100 seconds, and the standard 



ion 


<CF> 


166 seconds- 


In 


terms ot bins 


, the 


the 


standard devi at ion 


8- 


The standard 


devi « 


.imes 


the 


mean . 








5-1 


Number at Events 


per 


Interarr i val 


T i me 1 


li n 




Events 




Bin Events 


0 




67 








1 




36 




31 


0 


2 




19 




32 


0 


3 




15 




33 


1 


4 




12 




34 


0 


5 




7 




35 


0 


6 




4 




36 


0 


7 




3 




37 


0 


8 




6 




38 


0 


9 




1 




39 


0 


10 




3 




40 


0 


1 1 




4 




41 


1 


12 




2 




42 


0 


13 




1 




43 


0 


14 




3 




44 


0 


15 




2 




45 


0 


16 




0 




46 


0 


17 




0 




47 


0 


18 




0 




48 


0 


19 




0 




49 


1 


20 




1 




50 


0 


21 




0 




51 


0 


22 




0 




52 


0 


23 




2 




53 


1 


24 




1 




54 


0 


25 




1 




55 


0 


26 




1 




56 


0 


27 




2 




57 


0 


28 




1 




58 


0 


29 




1 




59 


0 


30 




0 
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Three different probability density functions were 
investigated to find an appropriate fit for Figure 5—3. 
They were 1) a hal f — gaussi an distribution, 2) an 
exponential distribution and 3) a J shaped distribution. 

The hal f -gaussi an probability density function isC33: 



p (t ) = 



~t/2t Q “ 



ysTftc 



(5-1 ) 



The general equations -for mean, mean square value and 
variance ( (j 2 ) can be used to solve -for the unkown 
constant, t Q : 






CO 

t p ( t ) dt 



■/ 



CO 



mean square value = J t p(t) dt 

0 

CO 

o / . . r> 



■I 



(J " = J (t - fl )"■ p(t) dt 
0 



= mean 



square value — 



(5-2) 



(5-3) 



(5-4) 



Substituting equation 5—1 into equations 5—2, 5-3 and 5—4 
leads to the -following relations: 



t V- , 



mean square value = t 



(5-5) 



G 






7 ( 



71 



t Q 2 ; (7= 0.756 jJi 
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This value -For t was used in equation 5-1, and the 
probability density -Function integrated over appropriate 
limits to get the number o-F events in each 20 second bin. 
The result is plotted against the experimental distribution 
in Figure 5—4. 

The second distribution (the exponential) belongs to 
the -Family o-F gamma distribution -Functions C153: 

i a -t/t Q 

t e . (5-6) 

< a + 1) 



t o a+ ip 



When Ql = 0, this becomes the exponential probability 
density -Function 



P 




-t/t, 



e 



(5-7) 



Again using equations 5—2, 5—3 and 5—4 leads to: 



2 

(X - t a ; mean square value = 2 fj, 

G 2 = /JL 2 ; <J = fl • 



(5-8) 



The exponential probability density -function was integrated 
over the bins, and the results are shown in Figure 5—5. 

Another demonstration o-f the -fit o-f the exponential 
probabilty distribution is shown in Figure 5—6. Taking the 
natural log o-f the -function should lead to a straight line 
when plotted against time or bin number . The straight line 



Number of Events 
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20 Second Bin (0=0- 19s, 1=20-39s, etc) 



Figure 5—4 HaH-gaussian distribution compared with 
experimental values. 



/ 



M lj m ber of Events 
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Figure 5-5 Exponential distribution compared with 
experimental values. 



C n u m fc> -e r of ovonts) 



-B^- 




Bin number 



Figure 5-6 Semi log plot o-f exponential distribution and 
experimental values against bin number. 
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m Figure 5-6 is a plot, o-f the natural log o-f the points 
calculated using the exponential probability density 
-function- The experimental points seem to curve rather 
than lie on a straight line. 

The last distribution <J shaped) is also a gamma 
distribution. The J shaped distributions are character i zed 
by Oi < 0 . I chose a -fairly common distribution with 
d = —0.5 . The probability density -function is: 

i ~ 1/2 -t/t D 

p ( t ) = t e , (5-9) 

and the key parameters are: 

' :,t o 

t Q = 2 ; mean square value = , (5-10) 

cr 2 - ; cr-V 2 ^ • 



This distribution is plotted against the experimental 
values in Figure 5-7. The natural log o-f both calculated 
and experimental paints are plotted against bin number in 
Figure 5-8. This distribution seems to -fit the 
experimental points best o-f all. The J shaped probabilty 
density -function goes to in-finity at zero, but it is 
i ntegrabi e. 

A Chi square goodness of -fit test was done on all 



three di str i butions. 



The results are summarized in Table 



M u m ber of Events 
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20 Second Bin (0=0- 19s, 1=20-39s, etc) 



Figure 5—7 J shaped distribution compared with 
ex peri mental val ues. 



eve nts) 
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Figure 5-8 Semi log plot o-f J shaped distribution 
experimental values against bin number. 



and 
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5—2. Also presented in Table 5-2 are the ratios of 



standard deviation to mean. 



Table 5-2 Comparison of Distribution Functions 



Chi square 



(J/fl 



Experimental 



1.66 



Hal f — gaussi an 



127. 17 



0. 76 



Exponenti al 



45.69 



1.00 



J shaped 



10.65 



1.41 



For a distribution to pass a goodness of fit test it 
must have a Chi square less than a prescribed limit. The 
limit for my test (9 degress of freedom, 01 = 0.005) was 
23.6C161. Only the J shaped distribution passed the Chi 
square test. It also has (j / jJi closest to the 
experimental values. In summary, the interarrival data 
reasonably fit a J shaped distribution given by: 



Since event detection rate depended on ambient noise 
level, interarrival time between events should also show 
environmental dependence. Table 5.3 gives average and 
standard deviation of the interarrival time for different 



p (t) 



i -1/2 -t/2 j! 



t e 



(5-11 ) 




ambient noise pressure levels. 
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Table 5—3 Background Noise Level Dependence o-f 

Interarrival Time 



Ambient Noise Interarrival Time 

rms pressure 

(20—80 Hz) mean standard deviation 



(Pa) 


(sec ) 


(sec) 


0.01-0.02 


67 


129 


0. 02-0.03 


190 


153 


0. 03-0.04 


147 


143 


over 0.04 


1S3 


318 



The tapes having a background noise level o-f 0-01 to 0-02 
Pa have a significantly shorter interarrival time than 
tapes in the other three pressure groups. As with 
detection rate, the in ter arr ival time does depend on 
ambient noise level . 

Spatial Analysis 

After removing nonlocatable events and events located 
outside a horizontal range of 20,000 meters, 164 events 
remained- These were grouped by horizontal range into 42 
annuli of equal area as shown in Figure 5—9- Each annulus 
is a 30 square km ring centered at the array origin. The 
first annulus ("0") went from 0 to 3090 meters, the second 
from 3090 to 4370 meters, and so on. 

Table 5—4 shows the number of events per annulus and 
Figure 5-10 shows this distribution graphically. 



I T l iouso nd 




Fi gur 



5-9 Scheme -for annuli each o-f area equal to 



Numfc>or of events 



100 



90- 
80- 
70- 
60- 
50- 
40- 
30- 
20 - 
10 - 
0 - 

0 10 20 30 40 

Range Annulus (30 square km each) 

Figure 5—10 Number a-f events per radius annulus. 
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Table 5—4 Number of Events per Annulus 



Annul us 


Events 


Annul us 


Events 


0 


91 


21 


0 


1 


19 




0 


n 


4 


97 


1 


7 




24 


0 


4 


2 


25 


0 


5 


4 


26 


1 


6 


9 


27 


0 


7 


3 


28 


0 


a 


1 


29 


0 


9 


5 


30 


0 


10 


*71 


31 


0 


i i 


"7 


-7- '■7 


0 


12 


4 


-j>- j» 


1 


13 


1 


34 


0 


14 


1 


35 


1 


15 


4 


36 


0 


16 


0 


37 


1 


17 


0 


38 


1 


18 


1 


39 


0 


19 


0 


40 


0 


20 


0 


41 


1 



The average number o-f events per annulus is 3.93 and the 
standard deviation is 14.15 events- Figure 5—10 shows that 
the number of events found is highly dependent on their 
range from the array. In the center annulus there were 
over 20 times the mean number o-f events. 

The dependence on range is not a surprise, since 
spreading (and possibly scattering and other losses) will 
reduce the strength of weak transients down to the ambient 
noise level. For this reason, the center annulus is 
probably the best indicator of actual event density. In 
this ring there were 91 events per 30 square kilometers per 
662 minutes of observation or approximately 0,3 events per 
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square kilometer per hour . 

The average number of events per annulus and the 
number of events per square kilometer per hour should 
depend on background noise level. The average number at 
events per annulus was -found for each ambient noise rms 
pressure range, and adjusted to reflect the number of 
events in a 662 minute period. The results are seen in 
Table 5-5. The number of events per square kilometer per 
hour for the center annulus are also shown in Table 5—5. 

Table 5—5 Average Number of Events per Annulus 
for 4 Ambient Noise Levels 



Ambient Noise 


Events 


Mi nutes 


Adjusted 


# Events 


rms pressure 


per 


of tape 


Events per 


per km^ 


(20-80 Hz) 


Annul us 


Examined 


Annulus 


per hr 


0.01-0.02 Pa 


2. 60 


287.5 


5.98 


0.452 


0.02-0.03 Pa 


0.29 


77 


2. 46 


0 . 260 


0-03—0.04 Pa 


0.45 


1 00 


2.99 


0. ISO 


over 0.04 Pa 


0 . 36 


157.5 


1 . 50 


0 . 076 


Ent i re 


F'opul ati on 


3.93 


662 


3. 93 


0.275 


The average 


number 


of events per 


annulus and 


the 


number of events 


per square kilometer 


per hour both reflect 


the effect of s i gn al —to- 


■noise rat i o o n 


the detect i 


on 



s c h em a . 

The entire population of events was investigated for 
angular dependence. Figures 5-11 and 5-12 show the number 



hJ in 



M u rn t> o r of ovonts 
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Figure 5-11 Number o-f events per 30° sector. Ancles dre 
measured from the northern leg of the hydrophone array. 
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000 




ISO 



Figure 5-12 Number o-f events per 30° sector- Radius 
gives the number o-f events, while angle indicates the 
sector measured -from the northern leg o-f the array, 
ring represents 10 events. 



Each 
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of events -found per 30° sector- In the polar diagram 
(Figure 5-12) the radius shows the number of events- The 
angles are measured from the northern leg of the array. 
Figure 5—13 is a polar plot showing the number of events 
per 10° sector. There was no predominant angular direction 
found . However, some preference can be seen for bearings 
of 330° and 190° from the northern leg of the array. 

Strength Analysis 

The mean hydrophone peak pressure magnitude for each 
event fell within a fairly narrow band of values. The mean 
peak pressures ranged from 1.32 to 0.16 Pa, with an average 
of 0.36 Pa and a standard deviation of 0.20 Pa. Figure 5—14 
shows the mean hydrophone peak pressure values for all 
events located between 100 m and 20,000 m plotted against 
range from the array origin. 

The different symbols shown in Figure 5—14 represent 
events during each of the four ambient pressure categories. 
The events w ith a higher mean hydrophone peak pressure have 
a tendency to occur during higher ambient noise levels - 
This can be seen in Table 5-6, where the maximum, minimum, 
average and standard deviation of the mean hydrophone peak 
pressure values are given for each of the four ambient 
noise levels. 
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Figure 5—13 Number o-f- events per 
gives the number of events, while 
sector measured from the northern 



10° sector. Radius 
angle indicates the 
leg o-f the array. 



- 100 - 




(Thousands) 

Horizontal Range from Array Origin (m) 

Q.O) -Q.Q2 Pa . Q Q.02-O.Q3fia. A 0.03 -O.OY Pa. X 0-04^ Pa. 



Figure 5-14 Mean hydrophone peak pressure measured for 
events between 100m and 20,000 m, plotted against 
horizontal range from the FRAM IV array origin. 
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Table 5-6 Mean Hydrophone Peak Pressure 
-for 4 Ambient Noise Levels 



Ambient Noise 


Mean 


Hydrophone 


Peak Pressure 




rms pressure 




(Pa) 






(20-80 Hz) 


max 


mi n 


average 


std dev 


0.01-0.02 Pa 


1.32 


0. 16 


0.31 


0. 17 


0.02-0.03 Pa 


0 . 56 


0.25 


0. 36 


0 . 09 


0.03—0.04 Pa 


0.85 


0.23 


0. 43 


0. 15 


over 0.04 Pa 


1 . 28 


0. 16 


0.49 


0.34 


Entire 










Popul ati on 


1.32 


0. 16 


0 . 36 


0 . 20 



Source strength (F Q ) was -found tor the events which 
had hydrophone locations between 300 m and 20,000 m from 
the event. 

The dipole strengths ranged from 33 kN to 4-9 MN , 
with an average of 431 kN and a standard deviation of 555 
kN- The distribution of strengths for the 151 events 
evaluated is shown in Figure 5-15 and in Table 5-7. 

Figure 5-16 shows the dipole strength for all events 
platted against horizontal range from the array origin. 
Again, it can be seen that the stronger events occur w hen 
the ambient pressure level is high . Table 5-8 gives the 
strength values for the different ambient noise levels. 



Table 5-7 Strength Distribution for a Population of Events 



F Q # of Events 



0 


to 


100 


kN 


19 


100 


to 


200 


kN 


29 


200 


to 


300 


kN 


33 


300 


to 


400 


kN 


28 


400 


to 


500 


kN 


8 


500 


to 


600 


kN 


5 


600 


to 


700 


kN 


4 


700 


to 


800 


kN 


5 


800 


to 


900 


kN 


X. 


900 


to 


1000 


kN 


5 


1000 


to 


1100 


kN 


3 


1 100 


to 


1200 


kN 


2 


1200 


to 


1300 


kN 


3 


1300 


to 


1400 


kN 


1 


1400 


to 


1500 


kN 


1 


1.5 


to 


3 


MN 


1 


over 




~T 


MN 


2 



Table 5-8 Dipole Strength 
Ambient Noise 



rms pressure 
(20-80 Hz) max 

0.01-0.02 Pa 1051 

0.02-0.03 Pa 2041 

0.03-0.04 Pa 1153 

over 0.04 Pa 4939 

Ent i re 

Population 4939 



versus Ambient Noise Levels 



Di pal e 


Strength (F_) 
(kN) 




mi n 


average 


std dev 


33 


259 


155 


64 


649 


634 


79 


643 


359 


59 


860 


1449 


33 


431 


555 



Nurriber of Events 
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Figure 5-15 Distribution o-f strength tor a population o+ 
noise events 
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(Thousands) 

Horizontal Range From Array Center (m) 



-j“ Q.QI-Q.oZPq. 0 OOZ'O-05 Pa. A 0.03'0-OVPa. X O-0H+- Pa. 



Figure 5-16 Dipole strength -For events between 300 m and 
20,000 m, plotted against horizontal range -From the array 
or i gi n . 
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Most of the events evaluated for strength occurred 
during the lowest ambient noise levels- Figure 5-17 shows 
the strength of events that occurred when the ambient noise 
was 0.01 to 0.02 Pa. The log of the dipole strength is 
plotted against the log of the horizontal range from the 
center of the array. The points scatter more so to the 
upper left rather than 1 ower right, because distance itself 
filters out weak events. A weak signal from far away would 
not reach the hydrophone array with enough amplitude to be 
distinguished from the background noise- And events 
located farther away would tend to be strong events. 
However, events located close to the array should have the 
entire range of source strength levels- This would produce 
a wedged shaped plot of weaker events close to the array. 
Indeed, Figure 5-17 shows a general scattering with perhaps 
a wedge of weaker events near the array origin. 

Nonetheless the trend shown in Figure 5—17 suggests 
that the ray average model used to estimate refractive- 
surface reflective spreading may need to be replaced with a 
more refined model. For example, horizontal ranges less 
than about 1000 m may include too small a loss, and 
therefore lead to too small a strength, because the 
reflective contr i but i ons may not be as large as imputed. 
Such a criticism is supported by the notion that for a 
given slope, £ , reflective rays and hence ray averaging 

occurs only beyond a critical horizontal range. 
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Figure 5-17 10 log of Dipole strength (dB re 1 N) versus 

10 log o-f horizontal range -from the array origin (dB re 1 m) 
-for events occurring during an ambient noise level of o ni 
to 0.02 Pa. 
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The -foregoing speculation suggests that the average 
dipole strength for the lowest- ambient noise case is best 
found from the events farther from the array, and is 

F Q a* 10 5 ' 5 N a* 320 kN , (5-12) 

with a much smaller standard deviation than in Table 5-8. 
Presumably corr espondi ng adjustments could be made for the 
higher ambient noise cases, but the FRAM IV data set 
contains too few events at higher ambient noise to plot as 
in Figure 5-17. 

The strength analysis is a somewhat ambivalent one 
because of spreading model uncertainty, and because data on 
ice slopes are not available. But the dipole picture of an. 
event likely has some validity, and at least rough 
estimates of its strength have been extracted from the 



data . 
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CHAPTER 6 

SUMMARY AND THOUGHTS 

Through the use of a detection program, visual 
confirmation and a location program, a population of 199 
Arctic noise transients was gathered. There are four major 
resul ts. 

First- more events are found when the ambient pressure 
is low, and more false alarms when the ambient pressure is 
high. The interarrival time and the average number of 
events per unit area also depend on ambient noise level. 
Since more events are found when the ambient noise is low. 
the interarrival time decreases, and the spatial density 
i ncr eases. 

Second, the interarrival times were fit to several 
possible probability distributions. The inter arrival time 
distribution best fits a J shaped gamma distribution. The 
mean inter arrival time is 100 seconds. 

Third, the number of events per unit area is highly 
dependent on range, since distance filters out weak 
transients- The event density in the annulus closest to 
the center of the array was 0.3 events per square kilometer 
per hour over all observations and 0-5 events per square 
kilometer per hour for quiet times. There is no 
predominant angular dependence to the spatial distribution 



of events 
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Last , the mean dipole strength for the observed events 
is 430 kN overall and 260 kN during low ambient noise 
levels- Stranger events occurred during high ambient noise 
levels- A refinement of the spreading loss model used to 
calculate these values may lead to values which are 
slightly higher. 

Analysis of Arctic acoustic events is far from 

complete- Several areas for improvement have been 
mentioned earlier in the thesis- The detection program 
needs to be made more robust to eliminate the event time 
error- A scheme for ignoring artifacts should be included. 
The location program wastes time looking in the wrong 
direction, although the bearing accuracy of the program is 
very good- The algorithm should be changed to quickly find 
the right bearing, and then search in a sector. 

The type of each event, whether it was a pop or a 
whine- was not recorded- Collecting this information and 
correlating it with interarrival time and range still needs 



to be done. 
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APPENDIX A 



User's Guide for the hdetect Program 



Figure A— 1: Flow Chart of the hdetect Program 



Source Code for the hdetect Program 
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USER'S GUIDE FOR THE hdetect PROGRAM 

The purpose o-f the hdetect program is to detect 
ambient noise transients amidst the background ambient 
noise recorded on a FRAM data tape- This is done bv 
comparing the short average o-f data paints to the long 
average o-f points on a single channel in order to -flag a 
possible detection, and then waiting until 50/- of the 
channels are -flagged to declare an actual detection. 

The input -for the hdetect program is a framread output 
•file without headers- A FRAM data tape is read into the 
file by the command 

-framread -head < RETURN > 

The program will ask for the input device (tape drive 
designation), the output file, and the number of data 
segments to skip and to read. Each segment represents 3.S 
seconds of data on 24 channels. The framread program reads 
a first segment which contains no data records, so you 
should specify skipping one more segment than you would 
normally calculate. For example, reading the entire first 
half of a 20 minute FRAM IV tape would require the response 
of 

1 160 < RETURN > 

to the question of "enter #skip, ^segment s : " . 

□nee this input file has been created the hdetect 
program can be used. The program is started with the 
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command 



hdetect 



RETURN > 



The program will ask for the FRAN tape number, the Julian 
date of the tape, and the start time o-f the tape in hours, 
minutes and seconds. The program will then ask you to 
select the channels you wish to use. In most cases the 
FRAN data tapes did not have ambient noise hydrophones tied 
into all channels, and the specific channel that a 
hydrophone was recorded on changed throughout the 
experiment. Which channels were in use and for which 
hydrophones can be found in the experiment logs. The 
program assumes that the channel number is equal to the 
hydrophone number, but allows you- to change thiis by 
inputting the channel number and the proper hydrophone 
number, or "0 M if the channel is not in use. For example, 
if channel 3 was not used, and channel 7 was used for 
hydrophone 21, the input would be 



The "0,0" ends the changes to the channel selection- You 
must now hit any key to continue the program. 

You will be asked to enter the input device (the input 
framread file), the number of skips and segments, and the 
name of the output file. The output file does not have to 
exist before the program is started. It will be created by 
the program. The number of skips and segments are those 



3,0 

7,21 

0,0 



< RETURN > 
< RETURN > 
< RETURN > 
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that you would calculate using 3.8 seconds per segment. 

For example, to proces the entire -first half of a 20 minute 

FRAM IV tape the number of skips and segments would be 

0,160 < RETURN > 

This is all of the input required by the user. The 

program proceeds from this point without user interaction. 

The output of the hdetect program is a file containing 

a list of detections in the following format: 

tapenumber Juliandate hour minute seconds 
0 eventnumber eventtime 

channel hydrophone timedelay amplitude 
channel hydrophone timedelay amplitude 
channel hydrophone timedelay amplitude 



channel hydrophone timedelay amplitude 
0 eventnumber eventtime 

channel hydrophone timedelay amplitude 

channel hydrophone timedelay amplitude 
0 eventnumber eventtime 

-1 

The "0" at the start of a line indicates a new event 
detection, and the — 1 11 at the start of a line indicates an 
end of file. Each channel that was flagged for a 
particular event is listed with its hydrophone number, 
timedelay from the earliest channel signal arrival, and its 
peak voltage amplitude. This outfile can be used as the 
input file for the location programs without modification. 
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Figure A-l Flow chart of the hdetect program. 
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Figure A-l Flow chart o-F the hdetect program. 
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Figure A-l 



Flow chart, o-f the hdetect program. 
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/* startdoc 
hdetect . c 

program to read whoi-segy format data tapes from the fram IV program, 
usage : 

hdetect 

program ia interactive. 



by Mary Townsend-Manning 

enddoc 

*/ 



♦include <stdio.h> 

♦include <math.h> 

♦define NCHAN 25 

♦define RECLN 950 /* number of samples per trace (4 bytes per sample)*/ 

♦define OBYTES 3800 /* number of bytes per record output */ 

♦define ZERO 0 

♦define LONGFILTLN 64 / * length of long average filter */ 

♦define FLN 64 
♦define MAXE VENTS 4 
♦define R£SET_DELAY 0.3 
♦define RATIO 2.38 
♦define SIG_DELAY 0.02 
♦define THRESHOLD 0.5 
♦define EVENT_DELAY 0.5 
♦define END -1 

double vconv(x,y,n) 

register float *x, *y; 

register int n; 



{ 

register double sum =» 0 . ; 



if (n > 0) 

{ 

do { 

sum +- *x++ * *y — ; 
) while { — n > 0) ; 

} 

return (sum) ; 



main ( ) 

{ 



Source cade -far the hdetect program 



float l_ave [NCHAN] , sh_ave, time of f lag [NCHAN] ; 
float ampof flag [NCHAN] , firsttime, timede lay [NCHAN] ; 
float chaneventtime [MAXEVENTS] [NCHAN] ; 

float chaneventamp [MAXE VENTS ] [NCHAN] , eventime [MAXEVENTS ] ; 

int flag [NCHAN]; 

int num_active__events ; 

int m, n; 

int event^f lag [MAXEVENTS] [NCHAN] , flag_sum, number_of_e vents ; 
int event^number [MAXEVENTS ] ; 
int nchan - 0; 

int chan [NCHAN] , k, j, i, tid, date, hour, min, sec, 1; 
int channel, data, toggle; 
char answer; 

float h [FLN] , longf ilt [LONGFILTLN] ; 
int nskip,nseg ; 

char oddobuf [NCHAN] [OBYTES ], evenobuf [NCHAN] [OBYTES ] ; 

int count, error; 

char fname [80] , iname [80] ; 

FILE *iptr, *ptr, *fp, *fopen(); 
float timeseries [NCHAN] [RECLN+2*FLN] ; 
float time * 0.; 

num__active_events - 0; 

for (i*l; i<25; i++) { 

chan [i] - (i) ; 
flag[i]- 0 ; 

for (m*l ;m<MAXE VENTS ;m++) { 

event_f lag [m] [i] » 0; 
chaneventtime [m] [i] = 0; 
chaneventamp [m] [i] = 0; 

} 

) 

/* Program initialization from keyboard */ 



fprintf (stderr, "Program Initialization\n" ) ; 
fprintf (stderr, "enter FRAM tape #\n") ; 
fscanf ( stdin, ”%d" , 5tid) ; 

fprintf (stderr, "enter Julian date\n") ; 
fscanf (stdin, "%d" , &date) ; 

fprintf (stderr , "enter time - HR, MN, SC\n" ) ; 
fscanf (stdin, "%d, %d, %d" , fihour, &min, &sec) ; 

fprintf (stderr, "default values for channels and phones\n" ); 
fprintf (stderr, "are channel # = phone #.\n"); 

fprintf (stderr, "enter channel, phone to change. \n") ; 
fprintf (stderr, "enter '0' for phone, to eliminate a channel. \n") 
fprintf (stderr, "enter '0,0' to quit.Xn"); 
fscanf (stdin, "%d, %d" , & j , 5k) ; 

while (j !- 0 55 j < 25) { 

chan[j] - k; 
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3 — 



f scanf (stdin, "%d,%d",&j,&k) ; 

) 

fprintf (stderr, "FRAM TAPE %d Julian Date: %d\n" , tid, date ) ; 

fprintf (stderr, " « Time: %d: %d: %d\n", hour, min, sec) ; 

for (i»l;i<13,* i++) 

fprintf (stderr, "CH %d PH %d CH %d PH %d\n", 

i, chan [i] , i+12 , chan [i+12] ) ; 



/* Check to make sure inputs are correct — Change if necessary */ 



f scanf (stdin, "%c" , & answer) ; 

fprintf (stderr, "Hit any key and RETURN, when ready."); 
f scanf (stdin, "%c" , & answer) ; 



for (i*l; i<NCHAN; i++) { 

if (chan [i] !» 0) nchan++; 

) 

fprintf (stderr, "enter input device: ") ; 

scanf ( " %s " , iname ) ; 

if((iptr - f open ( iname, "r")) — NULL) 

{ 

fprintf (stderr, "can' t open %s\n" , iname) ; 
exit (1) ; 

) 



fprintf (stderr, "enter #skip, tsegments : \n") ; 

fprintf (stderr, "values of 0 and 320 will read entire tape\n") ; 
scanf ("%d, %d", &nskip, &nseg) ; 

/* load bandwidth filter */ 

if ( (fp-fopen("PMf loat" , "r”) ) -= NULL) ( 

printf ("cannot open bandwidth filter file\n"); 
exit (0 ) ; 

) 

for (i*0 ; i<64; i++) 

fscanf (fp,"%f",&h[i]) ; 
fclose (fp) ; 



/* load averaging filter */ 



f or (i-0 ; i<LONGFILTLN; i++) 

longf ilt [i] » 1 . 0/ (float) LONGFILTLN; 



/* Open output file */ 
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fprintf (stderr, "enter out-file: ") ; 

scanf ( "%s" , fname) ; 

if ( (ptr - f open (fname, "w") ) — NULL) 

{ 

fprintf (stderr, ’’can' t open %s\n" , fname) ; 
exit (1) ; 

) 

fprintf (ptr, "%d %d %d %d %d\n", tid, date, hour, min, sec) ; 

time » time +3.8 * (nskip-1) ; 

/* enter first record */ 

if (nskip%2 == 1) { 

for ( j-1; j<NCHAN; j++) ( 

f read (ievenobuf [ j] [0] , sizeof ( float) , RECLN , iptr); 

) 

toggle - 1; 

) 

else { 

f or ( j-1 ; j<NCHAN; j++) { 

f read (Soddobuf [ j] [0] , sizeof (float) , RECLN, iptr); 

) 

toggle - 0; 

) 

/* ENTERING RECORD READING MODULE */ 
time - time + 0.504; 

/* fprintf (stderr, "using buffer size %d bytes\n" , sizeof (buf) ) ; */ 

for(i-l; i < nseg; i++) 

{ 

/* fprintf (stderr, "time - %f\n", time); 

fprintf (stderr, "processing record %d\n", nskip+i) ; */ 

/* read next record into appropriate buffer */ 

if (toggle ==» 1) { 

for (j-1; j<NCHAN; j++) ( 

fread (Soddobuf ( j] [0] , sizeof (float ) , RECLN, iptr); 

) 

toggle - 0 ; 

) 

else { 

for (j-1; j<NCHAN; j++) ( 

f read (Sevenobuf [ j] [0] , sizeof (float) , RECLN, iptr); 

) 

toggle - 1 ; 

) 

/* filter and square data */ 
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} 



if (toggle =» 1) { 

for ( j~l; j<NCHAN; j++) { 

if (chan ( j] !» 0) { 

sq_f ilt Uoddobuf ( j ] [0] , Sevenobuf ( j ] (0],h f 
^timeseries [ j] [0] ) ; 

l_ave[j] - vconv ( &longf ilt [ 0] , & time series [ j ] (2 * (FLN-1) ] , 
LONGFILTLN) ; 



} 

} 

else { 

for ( j=l ; j<NCHAN; j++) { 

if (chan [ j] !- 0) { 

sq_filt (Sevenobuf ( j] [0] , fioddobuf [ j] [0] , h, 

& times eries [ j] [0] ) ; 

l_ave(j] - vconv { &longf ilt [ 0 ] , & time series [ j ] [ 2 ^ (FLN-1 ) ] , 
LONGFILTLN) ; 

} 

} 



/* ENTERING EVENT DETECTION MODULE */ 

f or <k~0 ; k<R£CLN; k +« 4) { 

for (1=1; KNCHAN; 1++) { 

if (chan ( 1] !- 0) { 

l_ave[l] » ( 63 . 0*l_ave (1] + time series [1] (2* (FLN-1) +k] ) /64 . 0; 
sh_ave = (timeseries [1] (2* (FLN-1) +k] + 
timeseries [1] [2* (FLN-1) +k-l] + timeseries [1] [2* (FLN-1) +k-2] + 
timeseries [1] (2* (FLN-1) +k-3] ) /4 . 0; 

/* reset old flags */ 

if (flag [1] — 1 £& (time - time of f lag [ 1] ) > RESET_DELAY) 
flag [1] = 0; 



/* set flag if RATIO of signals is reached */ 



if ( (sh_ave/l_ave [1] ) >=RATIO) { 
if (flag (1] =»« 1) { 

if (sheave > ampof f lag [1] ) { 

timeof f lag [1] =* time; 
ampofflag(l] =» sh_ave; 

} 

} 

else { 

if (num_active_e vents =*= 0) { 

f lag [1] - 1; 
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timeof f lag [1] - time; 
ampofflag[l] * sh_ave; 

) 

else { 

for (m»l;m< (num_active_events + 1) ;m++) { 

if (event_flag [m] [1] -* 1) { 

if {{time - chaneventtime [m] [1] ) 

<- SIG_DELAY) { 

if (sh_ave > cha ne vent amp [m] [1]) { 
chaneventtime [m] [1] - time; 
chane vent amp [m] [1] * sh_ave; 

} 

} 

else { 

flag [ 1] - 1; 

timeof f lag [1] = time; 

ampofflagfl] - sh_ave; 

} 

} 

else { 



event_flag [m] [1] = 1; 
chaneventtime (m] [1] - time; 
chaneventamp [m] [1] * sh_ave; 
m * num active events; 



/* end of set flag module */ 

/* start new event module */ 

flag_sum = 0; 

for (1*1; KNCHAN; 1++) { 

if (flag[l] ** 1) flag_sum++; 

} 

if (( (float) flag_3um/ (float) nchan) >= THRESHOLD) { 
num_acti ve_events++ ; 
eventime [num_active_e vents] * time; 
number_of_events++; 

event_number [num^active^events ] * number_of_e vents ; 

for (1*1; KNCHAN; 1++) { 

event — flag [num_active — events] [1] * flagfl]; 
chaneventtime [num_active_events] [1] * timeofflag [1] ; 
chaneventamp [num_active_events] [1] * ampof flag [1] ; 
f lag [1] * 0; 
timeofflag[l] - 0; 
ampofflag[l] * 0; 
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} 



} 

/* end of new event module */ 

/* start of deactivate old event module */ 

if (num_active_events > 0 S& (time - eventime [1]) > EVENT_DELAY) { 
fprintf (ptr, "%d %d %f\n",ZERO, 

e vent_numbe r [ 1 ] , event ime [ 1 ] ) ; 

/* fprintf (stderr,"%d %f\n" , event_number [1] , eventime [1] ) ; */ 



/* find time delays by finding earliest channel event 

time, and subtracting that from the other channel 
times 

*/ 



firsttime = 10000.0; 

for (1=1 ; KNCHAN; 1++) { 

if{chan[l] !- 0 event_f lag [1 ] [1] !* 0 && 

chaneventtime [1] [1] < firsttime) 
firsttime =* chaneventtime [ 1 ] [1] ; 

} 

for (1-1; KNCHAN; 1++) { 

if (chan(l] !- 0 && event_ flag(l] [1] !** 0) { 

timedelay[l] - ( (chaneventtime [1] [1] ) - firsttime); 
fprintf (ptr, "%d %d %f %f \n", 1, chan [1] , 
timedelay [1] , chaneventamp [1] [ 1 ] ) ; 

} 

} 



/* print to file to indicate end of event */ 

for (1=1; l<num_active_events ; 1++) { 

for (m-1; m<NCHAN;m++) { 

event_flag [ 1] [m] * event_f lag [1+1] [m] ; 
chaneventtime [1] [m] - chaneventtime [1+1] [m] ; 
chaneventamp [ 1] [m] = chaneventamp [ 1+1] [m] ; 

) 

eventime [1] - eventime ( 1+1] ; 
e vent_n umber [1] - event_ number [1+1] ; 

1 

num_a c t i ve_e ven t s — ; 

} 

time +- 0.004*4.0; 

} 

} 

/* EXIT MODULE */ 
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/* fprintf (stderr, "processing record %d\n", nskip+nseg) ; */ 
if (toggle == 1) { 

for ( j=l; j<NCHAN; j++) ( 

if (chan ( j ] ! = 0) { 

for (k=0 ; k<OBYTES ; k++) 
oddobuf ( j] [k] = 0; 

sq_filt (fievenobuf [ j] [0] , fioddobuf [ j] [0] ,h, & time series ( j] [0] ) ; 

l_ave(j] = vconv (filongf ilt (0] , fitimeseries ( j ] (2* (FLN-1) ] , 
LONGFILTLN) ; 

} 

} 

) 

else { 

for ( j=l; j<NCHAN; j++) ( 

if (chan ( j] ! * 0) { 

for (k=0 ; k<OBYTES ; k++) 
evenobuf ( j ] [k] = 0; 

sq_f ilt (£ oddobuf [ j] (0] , fievenobuf ( j] [0] , h, fitimeseries ( j] [0] ) ; 

l_ave[j] - vconv (filongfilt [0] , & time series [ j] [2* (FLN-1) ] , 
LONGFILTLN) ; 

) 

) 

) 



for (k=0 ; k< (RECLN-2* (FLN-1) ) ;k+=4) ( 

for (1=1; KNCHAN ; 1++ ) { 

if (chan [ 1] != 0) { 

i_ave[l] = ( 63 . 0*l_ave [1] + time series (1] [2* (FLN-1) *k] ) /64 . 0 ; 
sh_ave = (timeseries (1] [2* (FLN-1) +k] + 
timeseries [1] (2* (FLN-1) +k-l] + timeseries [ 1] (2 * (FLN- 1) +k-2 ] + 
timeseries [1] (2* (FLN-1) +k-3] ) /4.0; 

/* reset old flags */ 

if (flag[l] == 1 && (time - timeof f lag [1] ) > RESET_DELAY) 
flag [1] = 0; 



/* set flag if RATIO of signals is reached */ 



if ( (sh_ave/l_ave [1] ) >=RATIO) ( 
if (flag ( 1] — 1) { 

if (sh_ave > ampof f lag [1] ) ( 

timeof flag [1] = time; 
ampof f lag [1] = sh_ave; 

) 

} 

else { 

if (num^active^e vents == 0) ( 

f lag [1] = 1; 
timeof flag [1] = time; 
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ampof f lag [1] * sheave; 

} 

else { 

for (m=l;m< (num_active_e vents + 1) ;m++} { 

if (event_f lag [m] [1] — 1} { 

if ((time - chaneventtime [m] [1}) 

<- SIG_DELAY) { 

if (sheave > chaneventamp [m] [1]) { 
chaneventtime [m] [1] - time; 
chaneventamp [m] [1] » sh_ave; 

} 

} 

else { 

flag [1] = 1; 

timeof f lag [1] = time; 

ampofflag[l] - sheave; 

} 

} 

else { 



event_f lag [m] [1] = 1; 
chaneventtime [m] [1] = time; 
chaneventamp [m] (1] = sh__ave; 
m » num active events; 



} 




/* end of set flag module */ 

/* start new event module */ 

f lag__sum = 0 ; 

for (1=1;1<NCHAN;1++) { 

if (flag[l] == 1) flag_sum++; 

} 

if (( (float) flag_sum/ (float) nchan) >= THRESHOLD) { 
num_a c t i ve_e vent s ++ ; 
eventime [num_active_events ] - time; 
n umb e r_o f_events++; 

event^number [ num_act i ve_e vent s ] = number_of_e vents ; 

for (1=1; KNCHAN; 1++) { 

event_f lag [num_active_events ] [1] =* flag[l]; 
chaneventtime [num_active_e vents ] [1] = timeof flagfl] 
chaneventamp (num_active_events ] (1] = ampof flag [1] ; 
f lag [1] = 0; 
timeof flag [1] = 0; 
ampofflag[l] =0; 
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/* end of new event module */ 



/* start of deactivate old event module */ 



if (num_ active_events > 0 && (time - eventime(l]) > EVENT_DELAY) { 
fprintf (ptr,"%d %d %f\n", ZERO, 

event_ number (1] , eventime (1] ) ; 

/* fprintf (stderr,"%d %f\n" , event_nunber [1] , eventime [1] ) ; */ 



/* find time delays by finding earliest channel event 

time, and subtracting that from the other channel 
times 

*/ 



firsttime - 10000.0; 

for (1-1; KNCHAN; 1++) ( 

if (chan [1] !- 0 && event_f lag [1] [1] !* 0 && 

chaneventtime ( 1] (1] < firsttime) 
firsttime * chaneventtime (1] [1] ; 

} 

for (1-1; KNCHAN; 1++) { 

if (chan [1] !- 0 && event_f lag [ 1] [1] !- 0) { 

timedelay [1] = ( (chaneventtime Cl] (1] ) - firsttime); 
fprintf (ptr, "%d %d %f %f \n " , 1, chan (1 ] , 
timedelay (1] , chaneventamp (1] [1] ) ; 

} 

} 



/* print to file to indicate end of event */ 

for (1=1; l<num_active__events ; 1++) ( 

for (m-1; m<NCHAN;m++) { 

event_flag(l] [m] - event_f lag [ 1+1 ] [m] ; 
chaneventtime [1] [m] - chaneventtime [1+1] [m] ; 
chaneventamp (1] (m] - chaneventamp [1+1 ] [m] ; 

) 

eventime [1] - event ime [1+1] ; 
event_number [1] - event_ number (1+1 ] ; 

1 

num_active_events — ; 

} 

time +» 0.004*4.0; 

) 

/* print out all events */ 
if (num_active_e vents >0) { 
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for (k-l;k< (num_active_events + l);k++) { 

fprintf (ptr,"%d %d %f\n", ZERO, 

event_numbe r [ k ] , event ime ( k ] ) ; 

/* fprintf (stderr,"%d %f\n", event_number (k] , eventime [k] ) ; */ 



/* find time delays by finding earliest channel event 

time, and subtracting that from the other channel 
times 

*/ 



firsttime - 10000.0; 

for (1-1; KNCHAN; 1++) { 

if(chan[l] !- 0 && event_f lag [k] [1] !- 0 && 

chaneventtime [k] (1] < firsttime) 
firsttime - chaneventtime [k] [ 1] ; 

) 

for (1-1; KNCHAN; 1++) { 

if(chan[l] != 0 && event_f lag [k] [ 1) 0) { 

timedelay [ 1] = ( (chaneventtime [k] [1] ) - firsttime) 
fprintf (ptr, "%d %d %f %f \n ,f , 1, chan [1] , 
timedelay [1] , chaneventamp [k] [1]) ; 

} 

} 

} 

} 

/* final summary to screen */ 



fprintf (ptr, "%d", END) ; 
fclose(ptr) ; 
fclose(iptr) ; 
exit(0) ; 






sq__f ilt (first, second, filter, output) 

/* filters and squares two data arrays */ 



float first [RECLN] , second [RECLN] , filter [FLN] ; 
float output [RECLN+2*FLN] ; 

{ 
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/* put zeros in first (FLN -1) places of output */ 
int j , i; 

float transition [2*FLN] ; 
float sum; 

f or (i*0 ; i<FLN-l; i++) 

output [i] * 0.0 ; 

/* put in first data */ 

f or (i-FLN-1; i<RECLN; i++) { 

sum « vconv < & f ilter [ 0 ] / & first [ i] / FLN) ; 
output [i] ~ sqrt (sum * sum) ; 

/* put in transition from first to second data */ 

for < j-0; j<FLN-l; j++) { 

transition [ j ] * first [RECLN- (FLN-1) + j ] ; 
transition [ j+FLN-1] - second[j]; 

} 

for (i-0;i<FLN-l;i++) { 

sum » vconv (&f ilter [ 0] , & transition [i+FLN-1] , FLN) 
output [ i+RECLN] - sqrt (sum * sum) ; 

) 

/* put in second data */ 

for <i-FLN-l;i<2*FLN;i++) { 

sum - vconv (&f ilter [0] , &second[i] , FLN) ; 
output [RECLN+i] ~ sqrt (sum * sum); 

) 



} 
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USER'S GUIDE FOR THE location? far locate? 
and finelacate PROGRAMS 

The purpose of the location programs is to find the 
spatial location of an event from the time delays between 
signal arrival at different hydrophones- The program 
assumes a test location and computes the slant range to the 
individual hydrophones- The slant ranges are plotted 
against the experimental time delays and a least squares 
fit is done. The test location with the best least squares 
fit is considered the location of the event. 

The input to the location programs is the output file 
of the detection program- Manual time delays may be 
substituted for the program generated time delays in this 
file, but this editting must be done before the location 
program is invoked. 

This location program is very user interactive. The 
user starts the location program with the command 
location < RETURN > 

The program asks for the input and output file names. It 
then reads the input file and asks whether the user would 
like to locate the first event. This allows the user to 
skip down to the event of interest. The program then 
allows the user to adjust which hydrophone time delays will 
be used in the location process. This is very hand'/ for 
removing questionable time delays, in order to get a better 
location solution. With the hydrophone channels chosen , 



the program proceeds with the actual location algorithm. 

In the finelcoate and far locate programs the user is asked 
to specify which quadrant or direction is to be searched. 

The location program tries test locations in a large 
grid, and when the "best" location is found, then searches 
a smaller grid around this "best " location. The location 
and far 1 ocate programs have 4 levels of grids and the 
f ine locate program has 3 levels. Intermediate answers are 
displayed for each level. 

The intermediate and final answers display the x and y 
coordinates of the best location, the standard deviation of 
the least squares fit (sigma), the group speed (which 
should be around 1440 m/s), and the y intercept of the time 
delay / slant range plot. After the final answer the user 
is asked whether or not he would like to remove outlying 
points. If this option is selected the program removes 
hydrophones with a deviation from the least squares fit of 
more than 3 times the standard deviation, and the group 
speed and the standard deviation are recalculated and 
di splayed . 

The user is then asked if he would like to locate the 
event with different hydrophones, and if so returns the 
user to the start of the channel selection process. The 
user may repeat this location scheme as many times as 
necessary for a particular event. Once the user is 
satisfied with the location answer , 
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the event with different phones, the program calculates the 
event strength based on the amplitudes in the input -File, 
the location o-f the event and spherical spreading losses. 
The location parameters and the strength are then written 
to an outfile. 

At this point the user is given the option to exit the 
program or locate the next event in the input file. The 
location process continues until the user exits or until 
the last event in a file is located. 
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/* location. c 

program to locate the spatial position of an event based on 
time delays taken from "detection”. 

Source strength is also computed. 

*/ 

# include <stdio.h> 

#include <math.h> 



^define PHONES 31 
#define LEVEL 4 
fdefine FINENESS 20 
#def ine DEPTH 91.0 

#def ine SENSITIVITY 0.0000000112202 



main () 



float amp [PHONES], phonex [PHONES ] , phoney [PHONES ] ; 

float timedelay [PHONES] , r [PHONES], bestrange [PHONES ] ; 

float sumtime, sumr, sumrsq, sumtimer, slope, yintr; 

float sigma, bestsigma, bestslope, bestyintr; 

float bestamp, N, source [PHONES] , sumsource; 

float xgs, ygs, xcntr, ycntr, a, b, bestx, besty; 

float level, gridsize, xfineness, yfineness, time, gpspeed; 

int i, j, tape, date, hour, min, sec, event, flag; 
int phone, nun, 1, n, m, bestflag; 
int phonef lag [PHONES] ; 

int eventselect, rerun, change, answer, bye, quadrant; 



char iname[80], oname[80]; 



FILE *ptr, open() , *locptr, *optr; 
/* PROGRAM INITIALIZATION */ 



bye - 0 ; 

fprintf (stderr, "input file =* \n") ; 
scanf("%s", iname) ; 

fprintf (stderr, "output file - \n M ); 
scanf("%s", oname) ; 



/* open files */ 

if ( (ptr * fopen (iname, "r")) NULL) { 
fprintf (stderr , "can' t open %s\n", iname); 
exit (1) ; 

) 
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if ( (optr = fopen(oname, "w") ) == NULL) { 
fprintf (stderr, "can' t open %s\n" , oname) ; 
exit (X) ; 

} 

if((locptr ■ fopen ("array_loc", "r") ) »- NULL) { 
fprintf (stderr, "can' t open array_loc file\n") ; 
exit (1) ; 

} 



/* read hydrophone locations into array */ 
for (i-1; i<P HONES ; i++) ( 

fscanf (locptr, "%d %f %f",&phone, &phonex(i], &phoney(i]); 

} 



/* read input file header */ 

fscanf (ptr, "%d %d %d %d %d", &tape, &date, ihour, &min, &sec) ; 
/* read event header */ 

fscanf (ptr, "%d ", &flag) ; 



while (bye !» -1) { 



eventselect * 0; 

while (eventselect !=* 1) { 

if (flag < 0) 
exit (0) ; 

for (i*l; i<P HONES ; i++) { 

phoneflag(i] = 0; 

} 

fscanf (ptr, "%d %f", ievent, &time) ; 

for(i=-l; i<P HONES ; i++) { 

fscanf (ptr, "%d" , &f lag) ; 

if (flag >0) { 

fscanf (ptr, "%d", &j); 
phoneflagtj] = 1 ; 

fscanf (ptr, "%f %f", &timedelay [ j ] , &amp(j]); 

} 

else { 

i - PHONES; 

) 

) 

fprintf (stderr, "event * %d, time = %f\n", event, time); 
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fprintf (stderr, "Do you wish to locate this event? (1 * yes)\n"); 
scanf("%d", fieventselect) ; 



/* channel selection */ 
rerun * 1 ; 

while (rerun =*=» 1) { 

change - 1; 

while (change — 1) { 

fprintf (stderr, " phone delay \n") ; 

for (i»l; i<P HONES ; i++) ( 

if (phoneflag (i] ! = 0) ( 

fprintf (stderr, " %d %f \n" , i, timedelay [i] ) ; 

} 

} 

fprintf (stderr, "Do you wish to change status? (1 = yes)\n”); 
scanf (”%d n , ^change) ; 
if (change « 1) { 

fprintf (stderr, "change status by typing phone#\n M ); 
fprintf (stderr, "type -1 to quit\n”) ; 
scanf ("%d" , & j) ; 

while (j !- -1) { 

if (phoneflag [ j ] !=» 0) { 

phoneflag [j] » 0; 

} 

else { 

phoneflag (j] * 1; 

} 

scanf ("%d", & j) ; 

} 

} 

) 



/* locate event */ 

num » 0 ; 
sumtime =» 0.0; 

for (i-1; i<PHONES; i++) ( 

if (phoneflag [i] !=» 0) ( 

num++; 

sumtime +=» timedelay ( i] ; 

} 

} 

fprintf (stderr, "select quadrant to search (1=NE, 2=*NW, 3=»SW, 4“SE\n") 
fprintf (stderr, "5-N, 6=S, 7*E, 8=W) \n") ; 
scanf ("%d", ^quadrant) ; 
if (quadrant =»= 1) ( 
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xcntr 


= 10000.0; 




ycntr 

} 

else if 


= 10000.0; 




(quadrant == 


2) { 


xcntr 


= -10000.0; 




ycntr 

) 

else if 


= 10000.0; 




(quadrant == 


3) { 


xcntr 


= -10000.0; 




ycntr 

} 

else if 


= -10000.0; 




(quadrant == 


4) { 


xcntr 


= 10000.0; 




ycntr 

i 


= -10000.0; 




i 

else if 


(quadrant == 


5) { 


xcntr 


= 0.0; 




ycntr 

\ 


= 10000.0; 




/ 

else if 


(quadrant == 


6) { 


xcntr 


= 0.0; 




ycntr 

\ 


= -10000.0; 




/ 

else if 


(quadrant == 


7) ( 


xcntr 


= 10000.0; 




ycntr 

} 

else if 


= 0.0; 




(quadrant == 


8) ( 


xcntr 


10000.0; 




ycntr 

} 


= 0.0; 




else ( 


xcntr 


= 0.0; 




ycntr 


= 0.0; 





} 



bestflag = 0; 

for (1=0; KLEVEL; 1++) { 

level = LEVEL-1-1; 
gridsize = pow ( 10 . 0 , level) ; 

for (m=0 ; m<FINENESS; m++) ( 

yf ineness = m - (FINENESS/2) ; 
ygs = ycntr + yfineness * gridsize; 

for (n=0 ; n<F INENESS ; n++) ( 

xf ineness = n - (FINENESS/2); 
xgs = xcntr + xfineness * gridsize; 
surar = 0.0; 
sumrsq = 0.0; 
sumtimer = 0.0; 
for (i-1; i<P HONES; i++) { 

if (phonef lag [i] !» 0) { 

a = xgs-phonex [i] ; 
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b =* ygs-phoney [i] ; 

r[i] * sqrt (pow (a, 2 .0) + pow(b,2.0) + pow (DEPTH, 2.0)) 

sumr +— r [i] ; 

sumrsq +- pow (r [i] , 2 . 0) ; 

sumtimer +- (r[i] * timedelay [i]); 



N = num; 

slope - ( (N * sumtimer) - (sumr * sumtime) ) / 
( (N * sumrsq) - pow (sumr, 2 . 0) ) ; 
yintr - (sumtime - (slope * sumr) ) /N; 



s igma - 0.0; 

for (i=l ; KPHONES; i++) ( 

if (phoneflag [i] != 0) ( 

sigma +- pow ( (timedelay [i] -yintr- (slope * r[i])),2.0) 

} 

} 

sigma * sqrt (sigma/N) ; 

if (bestflag — 0) ( 

bestx - xgs; 
be sty =* ygs; 
bestsigma - sigma; 
bestslope ■ slope; 
bestyintr = yintr; 
for (i-1 ; KPHONES; i++) ( 

if (phoneflag ! = 0) { 

bestrange[i] * r[i]; 

} 

} 

bestflag * 1; 



else ( 

if (sigma < bestsigma) { 
bestx = xgs; 
besty = ygs; 
bestsigma » sigma; 
bestslope * slope; 
bestyintr =» yintr; 
for (i-1; KPHONES; i++) ( 

if (phoneflag !- 0) { 

bestrangefi] =* r[i]; 

} 



) 



} 
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gpspeed =* 1 . 0/bestslope; 



fprintf (stderr, "bestx = %f, besty = %f, sigma =* %f\n", 
bestx, besty, bestsigma) ; 
fprintf (stderr, "group velocity * %f\n", gpspeed) ; 
fprintf (stderr, "y intercept = %f\n", bestyintr); 
xcntr - bestx; 
ycntr = besty; 



} 



fprintf (stderr, "Do you wish to remove outlying points? (l=*yes) \n") ; 
scanf ("%d", ^answer) ; 
if (answer == 1) ( 

nun = 0 ; 
sumr *0 . ; 
sumtime - 0 . ; 
sumrsq = 0 . ; 
sumtiner =■ 0.; 

for (i-1; i<PHONES; i++) { 

if (phonef lag [i] != 0) { 

, if (sqrt (pow ( (timedelay [i] -bestyintr- (bestslope *bestrange [i] ) ) , 
2.0)) < 2 . 5 *bestsigma) ( 
num++ ; 

sumtime +*» timedelay [i] ; 

sumr += bestrange ( i] ; 

sumrsq += bestrange ( i] *bestrange ( i] ; 

sumtimer +*■ bestrange [ i] *timedelay [i] ; 

} 

else { 

fprintf (stderr, "outlying Dhone # %d\n",i); 

} 

} 

} 

N = num; 

bestslope - ( (N*sumtimer) - (sumr*sumtime) ) / 

( (N*sumrsq) - (sumr*sumr) ) ; 
bestyintr = (sumtime- (bestslope*sumr) ) /N; 
gpspeed =* 1 . 0/bestslope; 

fprintf (stderr, "bestx = %f, besty » %f, sigma = %f\n", 
bestx, besty, bestsigma) ; 
fprintf (stderr, "group velocity = %f\n", gpspeed); 
fprintf (stderr, "y intercept = %f\n", bestyintr); 



fprintf (stderr, 

"Do you wish to relocate with different phones? (l=yes)\n"); 
scanf ("%d", irerun) ; 



} 
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/* finding source amplitude */ 

for (i*l; i<P HONES ; i++) { 

if (phonef lag [i] 1 “ 0) { 

sourcefi] » (amp {i] /SENSITIVITY) * bestrange [i] ; 
sumsource +*■ source [i]; 

} 

} 

bestamp = sumsource/N; 

fprintf (optr, "%d %d %f %f %f %f\n", event, date, time, 
bestx, besty, bestamp) ; 

fprintf (optr, "%f %f %f\n”, bestsigma, gpspeed, bestyintr) 

fprintf (stderr, "Do you wish to exit? (-1 = exit) \n") ; 
scanf("%d", &bye) ; 



f close (ptr) ; 
fclose (optr) ; 
fclose (locptr) ; 

exit ( 0 ) ; 

} 
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Table C-l 


Event 


Location Summary 






Name 




X 


y 


r 


R 


phi sigma 


s: gaa 


c Hydrophones 


ilaoe 


*. 














Removed 


event 


4 ! 


(iai 


(a) 


(ft) 


(a) 


(degrees) (sec) 


ta) 


(ii/sec) 



Na»c 15 the FRAh IV tape number followed by the time into the tape the event occurred, 

x and y are Cartesian coordinates, with the east leq of the array being the positive x axis. 

r is the horizontal range, and S is the slant range. 

phi is the bearing in degrees from the northern leq of the array, 

siqp.a * sec ) is the standard deviation of the time delays within an event. 

signs (a) is sigma (sec) tines the sound soeeo c. 

: :s the sourc ssee:. the inverse of the regression tit slope. 

t indicates that an event is considered non-locatable. 



4001,63 


16741 


-18733 


25310 


25311 


139 


0.0039 


5.62 


1425 5.17.22 


4001,126 


4413 


-5302 


6901 


6902 


140 


0.0034 


5.02 


1471 8,19,22 


4001,175 


-1041 


-9921 


9975 


9976 


136 


0.0014 


2.07 


1435 2,3 


4001,177 


-1380 


-10148 


10241 


10242 


188 


0.0058 


8.23 


1431 3 


4001.195 


5? g 


-2191 


2271 


2273 


165 


0.0034 


5.01 


1460 5.11.17 


4001,655 


at 3 


-2001 


2109 


2111 


lol 


0.0016 


2.30 


1444 22 


4001 ,315 


o417 


-9681 


11615 


11615 


146 


0.0019 


2.30 


1452 5,11,17 


4001.895 


-3352 


-10100 


10642 


10642 


193 


0.0025 


3.6? 


1470 none 


4001,106a 


3200 


-19330 


20136 


20136 


171 


0.0012 


1.77 


1478 17 


-001 . 1 131 


-9362 


9574 


13391 


13391 


316 


0.0027 


3.91 


1465 11.22 


J001.1156 


7 73 


-3101 


3196 


3197 


186 


0.0077 


11.27 


1457 5 


:.u=-4 


-1419 1 


!S3aO 


33204 


23205 


. Li. 


0.0043 


6.75 


!4 T 5 :: 


4 .'01 ,12 15 


2129 


-13591 


1901! 


19011 


1 7 4 


0.0074 


10. 94 


14' a none 


** ’ . , i I 


54 


-3*8 


358 


j / 0 


lee 


0.0134 


IS. 71 


1400 7 


4003,52 


-2321 


-9297 


9716 


9716 


197 


0.0039 


5.68 


1450 none 


4003.233* 


-216 


-30 


230 


243 


250 


0.0059 


8.14 


1368 18 


4003,254 


-1101 


-544 


1276 


1279 


240 


0.0093 


13.53 


1433 22 


4003,285 


265 


214 


341 


353 


51 


0.0108 


15.18 


1403 7 


4003,30! 


407 


145 


432 


442 


70 


0.0050 


7.11 


1435 13.18 


4003,419 


791 


4890 


4954 


4954 


9 


0.0026 


3.76 


1445 none 


4003.516 


142! 


-3001 


3320 


3322 


155 


0.0053 


7.64 


1434 none 


4003.634 


-6213 


4516 


7681 


7681 


306 


0.0033 


4.85 


1451 none 


4003.635 


2190 


-3001 


3715 


3716 


144 


0.0064 


9.10 


1420 7.12 


4003. ’03 


2890 


2230 


3650 


3652 


52 


0.0046 


6.54 


1407 12.13 


4-, 03. *95 


19995 


1802 


20076 


20076 


35 


0.0127 


IS. 26 


1443 i6,22.2- 


■*"v3. 1033 


.■it 


2-6 


J 


28" 


" r TS 

•J-J 


0.0036 


•J- . V 


1--2 P 
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Tacls C-l Event Location Suaaary 



Nase 


X 


v 


r 


A 


(Tape i. 










event 4) 


in) 


(a) 


(ft) 


(a) 


4005,219 


'312 


3890 


3902 


3904 


4005.241 


-103 


-178 


206 


226 


4005.395 


-342 


890 


953 


95a 


4005,803 


-88 


296 


309 


323 


4005.916 


-785 


616 


998 


1002 


4005 .997 


-15099 


-9810 


18006 


13006 


-005. i0! 2 


-1001 


-312 


2159 


2161 


4005,1097 


-4791 


-19503 


20033 


20083 


4005,1169 


1890 


-536 


1979 


1931 


4005,1173 


-49 


-140 


148 


175 


4005,1133 


-47 


-101 


111 


145 


4007,204 


323 


167 


364 


375 


4007.444 


791 


-108 


798 


304 


4007,582 


-205 


-74 


218 


237 


4007,794 


-197 


-73 


210 


230 


4007.320 


323 - 


485 


584 


591 


4007,938 


“-j’j 


-13 


jG 


100 



phi 


si gsia 


si Gina 


c 


Hydrophoi 

Removed 


(degrees) 


(sec) 


(in) 


(m/sec) 




355 


0.0107 


15.43 


1441 


none 


210 


0.0059 


8.39 


1430 


none 


339 


0.0239 


34.96 


1464 


7 


343 


0.0111 


16.09 


1450 


none 


308 


0.0037 


12.21 


1404 


none 


237 


0.0171 


24.51 


1433 


none 


248 


0.0327 


46 . 83 


’ 434 


■ one 


194 


0.0017 


2.50 


1447 


none 


107 


0.0039 


12.40 


1399 


none 


199 


0.0025 


3.57 


1409 


7 


205 


0.0010 


1.42 


1478 


4,9,10,1 


63 


0.0144 


21.21 


1469 


none 


93 


0.0053 


8.47 


1449 


none 


250 


0.0150 


20.93 


1398 


none 


250 


0.0156 


21.87 


1400 


none 


34 


0.0125 


17. 66 


1417 


- 


248 


0.0010 


1.54 


1488 


5 . c , 7 . . : 
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Tab! s C-l Event location 5usiaarv 



Naae 


x 


v J 
r 


r 


F; 


phi 


sigma 


sigsa 


c Bysr op hones 


(Taps 4, 
















Removed 


event 4! 


(ffl! 


(») 


(a) 


in) 


(degrees) 


(see) 


(si 


(a/sec) 


4009,100 


3390 


715 


3955 


3956 


80 


0.0041 


5.88 


1450 7,11,12,23,24 


4009,172 


-193 


313 


370 


332 


328 


0.0026 


3.74 


144? 23 


4009,173 


390 


427 


937 


992 


64 


0.0051 


7.40 


1464 15 


4009,193 


-2529 


920 


2735 


2737 


28? 


0.0046 


6.66 


14^8 none 


4009,194 


-537 


-1001 


1160 


1164 


210 


0.0032 


11.59 


1422 6,7,12,13,24 


4009,278 


3290 


-9023 


12253 


12253 


137 


0.0062 


8.85 


1432 none 




-i5 5 10 


•212" 


16031 


15051 


2 3 2 


0.0025 


3.52 


!42 7 ncr-3 


-■■09,300 


“C? 


150 


293 


303 


59 


0.0032 


4.57 


1423 13 


4009.320 


5553 


388? 


11046 


11047 


36 


0.0070 


10.21 


1455 none 


4009,334 


-407 


890 


979 


983 


335 


0.0118 


17,49 


1481 none 


1 009, 335 


-1575 


-3001 


3339 


3390 


200 


0.0051 


7.46 


1470 18 


4009,397 


3171 


131 


3174 


3175 


88 


0.0135 


19.49 


1441 none 


4009,410 


-175 


460 


492 


501 


339 


0.0071 


10. 10 


1414 none 


4009,447 


5939 


-3654 


7387 


7887 


118 


0.0022 


3.16 


1431 13 


4009,47a 


- 736 


488- 


333 


883 


304 


0.0049 


6.82 


1388 1,6.7,12 


4009,493 


773 


-2201 


2333 


2335 


161 


0.0032 


4,79 


1476 6,7,10,18 


4009.579 


-2551 


2612 


3651 


3652 


316 


0.0033 


4,74 


1457 none 


4009.600 


-13 


-35 


37 


100 


200 


0.002? . 


4.19 


1423 none 


4009,608 


23 


-1971 


1971 


1973 


179 


0.0043 


a. 17 


1423 none 


4009,320 


-1567 


2990 


3376 


3377 


332 


0.0047 


6,97 


1485 12 


4009,521 


;39 


-71 


693 


699 


96 


0.0057 


8.21 


1447 5 


400", 626* 


463 


574 


737 


743 


39 


0.0036 


13.15 


1534 19.20 


4009,5;" 


1265 


-2001 


2357 


236? 


148 


0.0275 


40.94 


U90 13.18 


-0 ) c , a” 7 


-m 


-1" T 


176 


19? 


219 


0.0031 


11.54 


1-36 none 


4309.725 


-745 


1390 


2032 


203.4 


538 


0.0090 


13.08 


1452 9.K- 


9009,785 


-1*2 


344 


335 


396 


■ j 0 - j 


0.0033 


5. 33 


[411 none 


4009.202 


~ ' j. ~ 


2390 


35=7 


3593 


j • 


0.0213 


31.46' 


;47S none 


4009,317 


739 


7902 


7934 


7934 


5 


0.007c 


to. SI 


1414 none 


4009,328 


390 


-a 76 


9 56 


971 


113 


0.0101 


14.32 


1424 none 


4009.373 


-160 


117 


198 


219 


306 


0.0153 


22.95 


1499 15.16.18,19,20 


4009,915* 


-23 


223 


229 


248 


354 


0.0122 


18.32 


1502 5,16 


4009,948 


-1639 


-3001 


341? 


3421 


209 


0.0070 


10.2? 


1473 13 


4009,94? 


3999 . 


-290 


9004 


9004 


92 


0.0036 


5.22 


1459 5 


4009,967 


-914 


-7510 


7565 


7556 


167 


0.0032 


4.62 


1452 none 


4009.932* 


105 


483 


499 


508 


12 


0.0143 


19.45 


1364 18,20,21 


4009,1041 


1390 


646 


1997 


2000 


71 


0.0107 


15.01 


1409 none 


4009.1150 


551 


153 


573 


501 


74 


0.0065 


9.23 


1416 18 
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Table C-i Even! Location Suflunsr 1 / 



Nase 


X 


V 


r 


R 


phi 


siosta 


51 Q ill a 


c Hydropnones 


(Tape i, 
















Resaved 


event #) 


ia) 


•;«) 


la) 


(mi 


(decrees) 


(sec) 


(a) 


(ffc/sec) 


4011,83** 


2890 


-699 


2973 


2975 


104 


0.0207 


27.41 


1326 none 


4011,168 


-922 


-6001 


607 1 


6072 


139 


0.0052 


7.61 


1462 none 


4011.169 


-6201 


5340 


3133 


3134 


311 


0.0060 


8.75 


1448 none 


4011,215 


-96 


-391 


403 


413 


194 


0.0042 


6.26 


1492 24 


4011,241 


2390 


-792 


2997 


2993 


105 


0.0118 


16.38 


1431 none 


4011.250 


-420 


890 


984 


989 


335 


0.0038 


5.6(1 


1491 7 


1,25? 


““ 1 

~ J -- i 




'4)3 


~r - r 


3 jo 


0.0035 


3.0’ 


‘ *rc r - - 

; 7 J j . ; • S 


401 1,278 


-584 


-1001 


1159 


lit; 


210 


0.0173 


24.27 


1401 7,11 


4011,203 


2190 


-301! 


3723 


3724 


144 


0.0032 


4.62 


1451 7 


4011,307 


3260 


5890 


6732 


6733 


29 


0.0060 


8.84 


1468 5,6,7 


4011,344 


1990 


766 


2132 


2134 


69 


0.0069 


10.29 


1431 15.18 


4011.362 


5573 


-10007 


11454 


11455 


151 


0.0097 


14.17 


1457 18 


4011,335 


-329 


-4001 


4086 


4087 


192 


0.0077 


11.22 


1463 none 


4011,510 


19870 


13790 


24186 


24187 


55 


0.0169 


23.52 


1393 5,6,7,19.20,21 


4011,532 


-561 


244 


612 


619 


294 


0.0068 


9.98 


1475 nans 


4011,541 


-258 


607 


660 


666 


337 


0.0084 


11.90 


1423 15,18 


4011,573* 


14471 


-16409 


21373 


2 137 5 


139 


0.0056 


8.31 


1578 6.15.16.19.20.2- 


4011.595 


358 


-32 


367 


37? 


103 


0.0102 


14.37 


i-OB 22 


4011,315 


3-212 


-4435 


'633 


i 333 


12a 


-0.0020 


2.32 


1433 : 3 , 13 


4011.316 


-530 


2890 


2958 


2959 


348 


0.0232 


33.55 


1447 none 


4011,663 


2990 


156a 


3375 


3377 


62 


0.0053 


7.81 


1479 18 


4011,634 


2:5 


256 


368 


380 


46 


0.0035 


4,87 


1393 22 


4011.693 


SO 


a 7 90 


379) 


-3791 


1 


0.0027 


3.33 


1452 21 


4011 , 7 C 9* 


3446 


• 334 


3a 1 2 


3613 


*3 


0.00 a 8 


9.41 


u'3 


4011.721 


-206 


-349 


7-7 

■J a. 


136 


220 


0.0132 


a a . . .' -j 


14;; none 


4011.811 


-396 


890 


974 


97? 


3;s 


0.0079 


11.75 


1-92 7 


4011,861 


136:3 


-15390 


20960 


20?60 


139 


0.0072 


10.1: 


i4ii 9.12.17.:: 


4011,863 


734 


-701 


1015 


1019 


134 


0.0112 


15.67 


1400 none 


4011,835 


1 90 


345 


1043 


1053 


71 


0.0053 


7.41 


140! 18.21.22 


4011.985 


6341 


-8212 


10375 


10376 


142 


0.0039 


5.61 


1423 none 


4011.992 


-1301 


1390 


2294 


2296 


325 


0.0238 


33.52 


1403 7,23,24 


4011,1046 


-655 


1390 


2000 


2002 


341 


0.0049 


7.26 


1487 7,12.13,23,24 


4011,1061 


-52 


107 


119 


151 


334 


0.0020 


2.70 


1333 3,4,5,6,20 


4011,1128 


- 636 


452 


780 


786 


305 


0.0047 


6.72 


1441 none 


4011,1137 


-3001 


-133 


3004 


3005 


267 


0.0097 


13.7? 


1416 23,24 
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Table C-l Event Location Summary 



Nate 


X 


V 


r 


R 


phi 


sigaa 


sigaa 


c Hydrophones 


(Tape #, 
















Removed 


event 0) 


<■) 


(a) 


(a) 


(a) 


(degrees) 


(sec) 


(a) 


(m/sec) 


4013,120 


-833 


-4970 


5039 


5040 


190 


0.0091 


13. 11 


1447 none 


4013,177 


370 


-1896 


1932 


1934 


169 


0.0035 


5.00 


1437 2,8,9,10 


4013,277 


-393 


366 


537 


545 


313 


0,0045 


6.38 


1418 8,9,14,16,22,23,24 


4013,317 


1326 


-12 


1326 


1329 


91 


0.0028 


3.99 


1445 none 


4013,351 


-54 


-99 


113 


146 


209 


0.0070 


9.72 


1396 none 


4013,354 


1534 


761 


1712 


1715 


64 


0.0016 


2.31 


1447 21 


4013,373 


216 


:83 


716 


722 


18 


0.0028 


4.11 


1450 none 


4013,331 


1248 


-2815 


3079 


3081 


156 


0.0051 


7.41 


1442 none 


4013,403 


2602 


-5450 


6039 


6040 


154 


0.0067 


9.56 


1420 none 


4013,464 


342 


2261 


2413 


2414 


20 


0.0083 


12.03 


1454 none 


4013,478 


-52 


662 


664 


671 


356 


0.0076 


10.56 


1388 none 


4013,519 


175 


492 


522 


530 


20 


0.0082 


12.05 


1472 14,16,17 


4013,661 


6390 


3658 


7801 


7801 


62 


0.0030 


4.34 


1463 none 


4013,664 


-1001 


1205 


1567 


1569 


320 


0.0114 


16.72 


1473 19,20,22 


4013,671 


19866 


-8441 


21585 


21585 


113 


0.0104 


14.72 


1416 none 


4013,694 


-675 


9290 


9314 


9315 


356 


0.0039 


5.62 


1437 none 


4013,723 


3100 


-1533 


3458 


3460 


116 


0.0059 


8.51 


1447 none 


4013.755* 


264 


168 


313 


326 


58 


0.0151 


23.00 


1520 24 


4013,776 


7989 


-5590 


9750 


9751 


125 


0.0052 


7.48 


1446 none 


4013,797 


-389 


7013 


7024 


7024 


357 


0.0052 


7.56 


1461 none 


4013,801 


-400 


-133 


422 


432 


252 


0.0078 


11.58 


1479 none 


4013,866 


38 


307 


309 


323 


7 


0.0062 


8.92 


1448 none 


4013,575 


49 


320 


324 


337 


9 


0.0031 


4.27 


1397 none 


*013,908 


9854 


17869 


20406 


20406 


29 


0.0066 


9.26 


1414 7,12 


4013,922* 


13 


26 


29 


97 


27 


0.0090 


13.61 


1514 8,9,10,11,12,18 


4013.950 


1165 


-2253 


2536 


2538 


153 


0.0037 


5.30 


1442 none 


4013,979 


7890 


920 


7943 


7944 


83 


0.0084 


11.99 


1432 none 


4013,1014 


-143 


2890 


2894 


2895 


357 


0.0091 


12.80 


1413 none 


4013,1018 


-137 


574 


590 


597 


347 


0.0169 


24.97 


1476 none 


4013,1107 


-701 


453 


835 


840 


303 


0.0081 


11.51 


1427 none 


4013,1124 


-9001 


8142 


12137 


12137 


312 


0.0064 


9.43 


1478 none 


4013,1166 


870 


-7001 


7055 


7055 


173 


0.0050 


7.39 


1468 18 


4013,1187 


-696 


693 


982 


987 


315 


0.0063 


9.34 


1472 17,18 


4013,1203 


7789 


-2315 


8126 


8126 


107 


0.0035 


5.11 


1480 none 


4015,232 


-101 


422 


434 


444 


347 


0.0026 


3.62 


1406 21,22 



4016 no events 
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Table C-l . Event Location Summary 



Name 


X 


y 


r 


R 


phi 


sigma 


sigma 


c Hydrophones 


(Tape t, 
















Removed 


event H) 


(at) 


(a) 


(a) 


(a) 


(degrees) 


(sec) 


(a) 


(fli/sec) 


4019,178 


-8416 


-18097 


19958 


19958 


205 


0.0038 


5.52 


1473 none 


4019,679* 


-17789 


10549 


20682 


20682 


301 


0.0023 


3.50 


1533 none 


4019,689 


-20109 


11680 


23255 


23255 


300 


0.0026 


3.70 


1428 none 


4019,776 


-21087 


12642 


24586 


24586 


301 


0.0037 


5.13 


1394 none 


4019,841* 


-3001 


-88 


3002 


3004 


268 


0.0061 


9.10 


1503 none 


4021,61 


-7356 


15898 


17608 


17608 


335 


0.0103 


14.68 


1430 none 


4021,154 


-6508 


18999 


20083 


20033 


341 


0.0062 


8.84 


1436 13 


4023 no 


events 
















4024 no 


events 
















4027,28 


-1001 


-371 


1068 


1072 


250 


0.0145 


20.81 


1434 24 


4029,666 


-4250 


-10010 


10875 


10875 


203 


0.0075 


11.20 


1489 3,6,11,13,29 


4029,897 


347 


-110 


364 


376 


108 


0.0244 


33.82 


1385 none 


4029,1059 


-3318 


-9110 


9695 


9696 


200 


0.0040 


5.71 


1439 29 


4031,823 


159 


- 373 


405 


416 


157 


0.0076 


11.32 


1493 none 


4033,490 


-2532 


4190 


4896 


4897 


329 


0.0047 


6.61 


1408 13,18,29,30 


4033,1041 


-21010 


767 


21024 


21024 


272 


0.0119 


17.80 


1493 9 , i 1,13,19,20,21 ,22,23.24 


4033,1043 


-10309 


1017 


10359 


10359 


276 


0.0054 


7.91 


1451 16,22,23,24 


4033,1053 


-285 


120 


309 


323 


293 


0.0042 


6.10 


1451 13 


4033,1126 


-21090 


1075 


21117 


21118 


273 


0.0020 


2.91 


1465 13 


4033,1173 


-20905 


805 


20920 


20921 


272 


0.0085 


12.89 


1515 13,29 


4040,26 


-20488 


19987 


28622 


28622 


314 


0.0125 


17.63 


1410 6 


4040,148 


-6001 


4252 


7355 


7355 


305 


0.0030 


4.32 


1426 10 


4040,385 


1462 


-4010 


4268 


4269 


160 


0.0054 


8.03 


1486 none 


4040,398 


-9001 


8110 


12116 


12116 


312 


0.0150 


21.35 


1427 13 


4040,714 


-5432 


-14089 


15100 


15100 


201 


0.0056 


8.38 


1497 none 


4040,730 


6468 


-17109 


18291 


18291 


159 


0.0126 


17.76 


1411 6 


4040,871 


-3424 


-8910 


9545 


9546 


201 


0.0022 


3,20 


1446 none 


4040,1114 


-20973 


16769 


26853 


26853 


309 


0.0061 


9.03 


1469 11 


4040,1157 


-7408 


7391 


10823 


10824 


317 


0.0075 


11.26 


1493 none 


4047,984 


-457 


200 


499 


507 


294 


0.0230 


31.93 


1387 none 


4049,1064 


-2320 


-8110 


8435 


8436 


196 


0.0068 


9.75 


1428 none 


4049,1066 


-1889 


-8001 


8221 


8221 


193 


0.0150 


21.01 


1401 9,11,13 
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Table C-l Event Location Suaaary 



Nase 




X 


y 


r 


R 


phi 


sigma 


sigma 


c 


Hydrophones 


(Tape #, 
event f) 




(M) 


(a) 


(a) 


(a! 


(degrees) 


(sec) 


(a) 


(m/sec) 


Removed 


4051 


no 


events 


















4053,3 




-16878 


18965 


25388 


25388 


318 


0.0021 


3.11 


1461 


7 


4055 


no 


events 


















4057 


no 


events 


















4059 


no 


events 


















4061,571 




-8393 


-20923 


22544 


22544 


202 


0.0127 


18.29 


1445 


21,22,23 


4061,379 




-3355 


-2455 


4157 


4158 


234 


0.0027 


3.91 


1445 


5,6,23 


4063,296 




■10503 


17985 


20827 


20827 


330 


0.0090 


12.91 


1441 


none 


4063,824 




■16793 


14787 


22375 


22376 


311 


0.0043 


6.29 


1455 


none 


4063,864 




-2001 


1346 


2412 


2413 


304 


0.0043 


6.18 


1424 


18 


4063,933 




-9077 


7976 


12083 


12084 


311 


0.0111 


15.60 


1408 


none 


4063,936 




-1130 


912 


1452 


1455 


309 


0.0037 


5.42 


1450 


18 


4065 


no 


events 


















4067 


no events 


















2001,543 




-3455 


-5112 


6170 


6171 


214 


0.0026 


3.77 


1452 


o 


2009,139 




-1997 


-20937 


21082 


21082 


135 


0.0036 


5.40 


1487 


none 


2009,453 




374 


-5201 


5274 


5275 


170 


0.0070 


9.30 


1402 


18 


2009,455 




1441 


-2001 


2466 


2468 


144 


0.0045 


6.55 


1445 


12 


2009,925 




-105 


-o73 


681 


687 


189 


0.0193 


27.43 


1422 


none 


2023,74 




15430 


-16460 


22561 


22562 


137 


0.0104 


15.13 


1458 


7,11 


3001,11 




-3921 


75 


3922 


3923 


271 


0.0125 


17.92 


1434 


none 


3001,16 




-20103 


-9471 


22222 


22222 


245 


0.0137 


19.92 


1450 


none 


3001,301 




-3001 


698 


3081 


3083 


283 


0.0149 


21.13 


1418 


4,5,6,7,20,21 ,22,2 


3047,169 




-19109 


-2181 


19233 


19233 


263 


0.0073 


10.56 


1453 


23 


3047,751 




890 


-538 


1040 


1044 


121 


0.0134 


19.52 


1460 


none 












N = 




199 


199 


199 














mean = 




0.0077 


11.09 


1444 














std dev - 




0.0055 


7.82 


34 
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T a g 1 e C-2 Tape Sugary 



TiOS i 


Minutes 


Julian Detections Art i i acts 


Events 


Multiple 


False 


P ras 




Examined 


Date 








Responses 


Alans 


(20-80 H: Pa) 


4001 


20 


S6 


35 


0 


13 


5 


17 


Net Aval 1 


2001 


20 


38 


8 


0 


1 


0 


c 


Not Avail 


2009 


20 


89 


15 


0 


4 


1 


10 


0.022 


3001 


so 


90 


5 


0 


7 


0 


n 

s. 


0.013 


4003 


20 


91 


15 


0 


13 


1 


1 


0.044 


4005 


20 


91 


24 


5 


11 


7 


6 


0.035 


4007 


17 


91 


11 


i) 


6 


I 


4 


0.022 


4)09 






55 


7 


3? 


b 


T 


0.1 10 


401! 




56 


47 


7 


35 


1 


4 


0 . 0 1 o 


4013 


20 


93 


48 


i 

i 


34 


5 


a 


0.014 


2023 


20 


98 


1 


0 


1 


0 


0 


0.037 


4015 


20 


99 


14 


8 


1 


0 


5 


0.040 


3047 


17.5 


103 


2 


0 


7 


0 


0 


0.010 


401a 


20 


105 


3 


3 


0 


0 


0 


0.017 


4019 


20 


105 


14 


3 


5 


1 


5 


0.016 


4021 


20 


109 


6 


0 


7 

i. 


7 

L 


2 


0.013 


4023 


20 


109 


22 


20 


0 


0 


4 


0.011 


4024 


20 


109 


14 


li 


0 


Ci 


3 


0.011 


4027 


20 


L0 


wC 


■-•C 


1 


ij 


A 


0.012 


4029 


20 


no 


11 


8 


3 


0 


0 


0.012 


4031 


20 


110 


11 


10 


1 


0 


0 


0.012 


4033 


20 


110 


19 


12 


6 


0 


1 


0.012 


4040 


20 


111 


•v* 


to 


9 


1 


3 


0.034 


404' 


20 


4 I 1 

ill 


6 


0 


1 


0 


c 


0.114 


4049 


20 


Li 


j 


0 








0. 142 


4051 


20 


ill 


0 


0 


0 


0 


0 


0.140 


40:3 


20 


112 


’■•I 
i. i 


0 


1 


7 


10 


0.080 


4)55 


r.s 


n: 


0 


0 


0 


o 


0 


0. '/Pi! 


4057 


20 


112 


9 


1 


0 


0 


8 


0.053 


4059 


20 


112 


1 


0 


0 


0 


1 


0.065 


4061 


20 


112 


6 


0 


7 

i. 


0 


4 


0.034 


4063 


20 


112 


9 


0 


J 


3 


1 


0.023 


4065 


20 


112 


0 


7 

L 


0 


0 


\ 


0.027 


4067 


20 


112 


6 


6 


0 


0 


0 


0.013 


Total 


662 




499 


139 


199 


36 


125 






Datscticr. 


is declares 


dhsn st least 


507. 


the hvdroohor.es have 


t; sages. 






Artifacts 


are san-nade 


noises -'shots, air gun blasts) 










i ■ S't s are 


decscticns 


strong enough 


to anal 


y:s. 









Multiple Responses are multiple detecticns of an already counted event. 
False alarms are detections too weak to analyze. 
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Tibi? C-3 Event Interarrival Time Summary 



Nase 


Event Time 


Inter arrival 


Time 


Han cental Rang 


(Tase 1, avsnt 


f) (sec) 


(sac) (20s 


bin) 


(a i 


4001,63 


63 


63 


3 


25310 


4001,126 


125 


63 


•y 

•J 


6901 


4001,175 


175 


49 


2 


9975 


4001,177 


177 


2 


0 


10241 


4001,175 


195 


18 


0 


2271 


4001,658 


653 


463 


23 


2109 


*■'01.318 


318 


160 


3 


11615 


i:01.595 


£95 


7" 


■-> 


10642 


4001 , iOit 


1066 


171 


8 


20136 


4001,1131 


1131 


65 


3 


13391 


4001,1156 


1:56 


25 


1 


3196 


4001 .ll c 4 


1194 


38 


1 


23204 


4001,1215 


1215 


21 


1 


19011 


4003,21 


21 


21 


1 


358 


4003.52 


52 


31 


1 


9716 


4003,293 


233 


231 


11 


230 


4003,234 


284 


1 


0 


1276 


4003.235 


235 


\ 


0 


34 1 


4003.30! 


301 


16 


0 


432 


4003.419 


41? 


113 


r 

•J 


4954 


4003.51s 


516 


97 


4 


3320 


4003.634 


634 


118 


5 


7681 


if 03.-35 


635 


1 


0 


3715 




'08 


■7'* 


*• 


36:0 


4003.395 


995 


237 


14 


I'M} 7 a 


4003.1053 


1053 


58 


7 


272 


4005.219 


219 


21 Q 


10 


3902 


4005,24! 


241 


22 


1 


206 


4005.395 


395 


154 


7 


953 


4005.303 


308 


413 


20 


309 


4005, 9 ! 0 


916 


108 


j 


998 


4005.997 


997 


31 


4 


13006 


4005.1012 


1012 


15 


0 


2159 


4005,1097 


1097 


35 


4 


20083 


4005.1169 


1169 


72 


3 


1979 


4005.1173 


1173 


4 


0 


143 


J' '5.1163 


1 183 


10 


0 


tr. 




L'j 4 


204 


iC 


3-jV 


**:Q7,444 


444 


240 


12 


798 


4007,552 


532 


133 


6 


218 


4007,794 


"94 


212 


10 


210 


, • • -» - T • 




2 b 


: 


5S4 




9 5b 


led 


3 


-.0 
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Tab I a C-3 Event Interarrival Tune Suasary 



Naae 


Event Ti®e 


Interarrival Ti&e 


Horizontal Rani 


(Tape it, event * 


) (see) 


1 sect 


<20s bin) 


(mi 


4009,100 


100 


100 


5 


3955 


4009,172 


172 


-*n 


3 


370 


4009,173 


178 


6 


0 


987 


4009,193 


193 


15 


0 


2735 


4009,194 


194 


1 


0 


1160 


4009,278 


27 B 


84 


4 


12253 


4009.292 


232 


14 


0 


16051 


4 > 09,200 


300 


a 


0 


233 


4009.350 


380 


60 


4 


11046 


4009,334 


334 


4 


0 


979 


4009,335 


335 


1 


0 


3339 


4009.397 


397 


12 


0 


3174 


4009,410 


410 


13 


0 


492 


4009,447 


447 


37 


1 


7387 


4009,476 


476 


29 


1 


383 


- 4009,493 


493 


17 


0 


2333 


4009,579 


579 


86 


4 


3651 


4009,600 


600 


21 


1 


37 


4009.608 


603 


8 


0 


1971 


4009, o20 


a 20 


12 


0 


3376 


4009.621 


621 


i 


0 


693 


4009.626 


626 


5 


0 


737 


4009,667 


667 


41 


n 

i. 


2367 


4009. a 77 


67' 


10 


0 


176 


4009,728 


723 


31 


- 


2032 


4009,78s 


736 


58 


i. 


385 


4005.302 


802 


16 


0 


3597 


4009,317 


81' 


15' 


0 


7934 


4009.825 


823 


11 


0 


96o 


4009,378 


373 


50 


•1 

i. 


133 


4009.916 


916 


38 


1 


229 


4009,948 


943 


32 


1 


3419 


4009,949 


°49 


t 


A 


9004 


4009,967 


967 


18 


0 


7565 


4009.932 


982 


15 


0 


499 


4009,1041 


1041 


59 


2 


1997 


4009,1160 


1160 


119 


ti 


573 
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r iDle C-3 Event Intararrivai Tiae Suasiary 



Nase 


Event Ties 


inter arrival 


Tiiie 


Horizontal Ran! 


(Tape 5 , event si 


!isc! 


i sec ) {205 


bin) 




4011,33 


83 


83 


4 


2973 


4011,168 


163 


35 


4 


c-071 


4011,16? 


169 


1 


0 


8183 


4011,215 


215 


46 


0 

L 


403 


4011.241 


241 


26 


1 


2997 


4011.250 


250 


? 


0 


904 


4011,25? 


259 


9 


0 


7908 


-011.278 


273 


t? 




.159 


4011,283 


283 


5 


0 


7723 


4011,307 


307 


24 


1 


6732 


4011,344 


344 


37 


1 


2132 


4011,362 


362 


18 


0 


11454 


4011.335 


385 


23 


1 


4086 


4011,510 


510 


125 


6 


24186 


4011,532 


532 


Li. 


1 


612 


4011,541 


541 


9 


0 


660 


4011,573 


573 


32 


1 


21878 


4011, 5?5 


595 


22 


1 


36 / 


4011.ol5 


2:5 


20 


< 

i 


' 6 o-j 


4011.616 


616 


i 


0 


2958 


4011,623 


6 63 


47 


0 

L 


3375 


4011,664 




1 


0 


368 


40!l,t53 


6 3 3 


25 


1 


6790 


1011. '0 5 


'Os 


12 


0 


r ,r ' ^ 


4011,211 


511 


90 


4 


974 


4011.861 


82! 


50 


n 


20960 


4011,363 


363 


* 

L 


0 


1015 


4011.305 


835 


22 


1 


1048 


4011,985 


985 


100 


5 


10375 


4011,992 


992 


7 


0 


2294 


4011,1046 


1046 


54 


0 

L 


2000 


4011,1061 


1061 


15 


0 


119 


4011.1123 


1128 


67 


3 


730 


4011,1137 


1187 


59 


0 

L 


3004 
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Tafcie C-3 Event Interarnval lias Suausarv 



Naas 


Event Tiffie 


Interarrival 


Tiffie 


Horizontal Rani 


(Taps #. event t! 


(sec) 


(see! (20s 


bin) 


■.a) 


4013.120 


120 


120 


6 


5039 


4013.177 


177 


57 


2 


1932 


4013,277 


277 


100 


5 


537 


4013,317 


317 


40 


2 


1326 


4013,351 


351 


34 


1 


113 


4013.354 


354 


3 


0 


1712 


4013.373 


373 


19 


0 


716 


- 13.3a. 


:ai 


s 


0 


307= 


4013.403 


403 


22 


1 


6'v39 


4013,464 


464 


61 


3 


2413 


4013,473 


478 


14 


0 


664 


4013,518 


j13 


40 


2 


522 


4013,661 


661 


143 


7 


7801 


4013,664 


6s4 


7 


0 


1567 


4013,671 


671 


7 


0 


21585 


4013,694 


694 


23 


1 


9314 


4013,723 


723 


29 


1 


3458 


4013,753 


735 


t n 
Ji 


1 


313 


40 13. ’76 


77c 


2! 


1 


9750 


4013.797 


797 


21 


i 


7024 


4013.801 


801 


4 


o 


422 


4013,36s 


866 


65 


■j 


309 


4013.375 


373 


9 


0 


324 


4013.903 


-08 


’7 


1 


20406 


4013.422 


422 


14 


o 


:? 


4013.350 


950 


23 


1 


cJvj 


4013. 9'3 


979 


29 


1 


7 5 43 


J 0.3. ,01- 


1014 


CJ 


1 


2394 


4013,1018 


1018 


4 


0 


590 


4013,1107 


1107 


89 


4 


835 


4013,1124 


1124 


17 


0 


12137 


4013,1166 


1166 


42 


2 


7055 


4013,1187 


1107 


21 


i 


982 


4013.1203 


1203 


16 


0 


8126 


4015,232 


232 


232 


11 


434 



401b 



ncne 
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Table C-3 Event Interarrival Time Summary 



Naie 


Event Tiee 


Inter arrival 


Time 


Horizontal Range 


ITaoe I, event # 


) (sec) 


(seci (20s hin) 


(ffl) 


4019,178 


178 


17E 


8 


18958 


4019,679 


679 


501 


25 


20632 


4019,669 


68? 


10 


0 


23255 


4019,776 


77a 


87 


4 


24536 


4019,841 


341 


65 


■j 


3002 


4021,81 


81 


81 


4 


1760S 


-'■1.154 


134 


73 


■j 


20083 



4023 none 

4024 none 



4027,23 


23 


28 


i 


1068 


4029,666 


666 


666 


33 


10875 


4029,897 


897 


231 


11 


364 


4029,1059 


1059 


162 


8 


9695 


4031.823 


377 


323 


*1 


405 


4033,480 


490 


490 


24 


4396 


4033.1041 


1041 


CC 1 
JJ i 


21 


21024 


4033,10*3 


1048 


7 


0 


10359 


4033.1053 


1053 


c 

J 


0 


309 


*033. 112: 


112a 


i 




? 1 11 7 


*033,1173 


1173 


47 


- 

L 


20920 


40*0,36 


j.6 


26 


l 


23622 


4040,148 


146 


122 


u 


7355 


4040,335 


335 


237 


11 


4268 


4040,398 


393 


13 


0 


12116 


4040,714 


714 


316 


15 


15100 


4040,730 


730 


16 


0 


18291 


4040,371 , 


371 


141 


7 


9545 


4040,1114’ 


1114 


243 


12 


26853 


4040.1157 


1157 


43 


2 


10323 


4047.934 


934 


984 


49 


499 


*0*9.1064 


1 i;LS 


1064 


CT 
w j 


8435 


*■ ** , ;/aa 


1 0 0 C 


• 


0 


8221 
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Table C-3 Event Interarrival Tise Summary 



Na*e 


Event Time 


Interarri val 


Ti me 


Horizontal Ran-; 


^Taoe 1, event 


4) (sec) 


(sec) (20 


5 bin) 


(mi 


4051 


none 








4053,3 


3 


3 


0 


25388 


4055 


none 








4057 


none 








405- 


none 








-Oil .571 


57! 


571 


20 


22544 


4061,879 


379 


308 


15 


4157 


4063,296 


29* 


2?6 


14 


20027 


4063,824 


824 


523 


26 


22375 


4063,364 


364 


40 


L 


2412 


4063,933 


933 


69 


3 


12033 


4063,936 


936 


3 


0 


1452 


4065 


none 








4067 


none 








2001,543 


543 


543 


27 


6170 


2009,189 


189 


1 rjc 
iu ' 


9 


21082 


2 C 0 9 .453 


453 


26“ 


« 7 


3274 


2:09,45: 


452 




•i 


2-66 


2009, -"15 


925 


470 


Z-j 


631 


2023.74 


7 4 


74 


j 


2256 1 


3001,11 


11 


11 


0 


3922 


3001,16 


16 


5 


0 


n'tnnn 

i-i-LLL. 


3001,301 


301 


285 


14 


3081 


3047,169 


169 


169 


3 


19233 


3047,731 


75! 


582 


29 


1040 




N = 


199 








'=*- = 


100 








at i oev - 


166 
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Taole C-4 Event Strength Summary 

Name r Po si qma Fc sigma f 

(Tape s, event #) 1st) (Pa) (Pa) (N) (N) (Hz) 

Naas is the FRAil IV tape number followed by the time -into the tape the event occurred, 
r is the horizontal range. 

Po is the average peak hydrophone pressure tor an event. 

sigma (Pa) is the standard deviation of the peak hydrophone pressure within an event. 

Fo is the average peak dipole strength calculated for an event. 

sigma iH) is the standard deviation of the dipole strength within an event. 

f is the frequency, explicitly, the axis crossing rate. 



4001 , 126 


6901 


0.487 


0.092 


707540 


134109 


53 


4001,175 


9975 


0.600 


0.108 


966109 


175081 


4e 


4001,177 


10241 


0.545 


0.194 


1002782 


357437 


42 


4001.195 


2271 


0.558 


0.117 


896709 


185120 


42 


4001.653 


2109 


0.665 


0.097 


1340877 


208247 


32 


4001,318 


11615 


1.117 


0.213 


1418008 


270370 


61 


4001,395 


10642 


0.524 


0.085 


1094608 


177188 


37 


4001,1131 


13391 


0.482 


0.085 


779610 


136750 


4S 


4001,1156 


3196 


0.495 


0.069 


1222273 


166340 


30 


4001,1215 


19011 


0. 3 27 


0.101 


1229059 


198513 


40 


4003,21 


j 3 


0.233 


0.103 


58368 


13452 


71 


4003,52 


3716 


0.324 


0.063 


430173 


83712 


58 


4003,234 


1276 


0.291 


0.131 


187664 


30580 


64 


4003,285 


34! 


0.362 


0. 136 


61870 


14917 


BO 


4003.301 


432 


0.274 


0.057 


62054 


15905 


74 


4003.4:: 


4?f4 


3. 436 


0.450 


622783 


640604 


53 


4003.516 


3320 


0.246 


0.046 


342405 


63384 


5-j 


4 v 0 ■ 3^4 


763! 


0.163 


0.048 


235652 


69027 


CT 
O v' 


4003.335 


3715 


0.242 


0.044 


368396 


66290 


49 


4003,703 


3650 


0.225 


0.063 


264688 


7373! 


64 


4005,219 


3902 


0.305 


0.271 


358094 


314253 


64 


4005.395 


953 


0.312 


0.170 


158324 


70515 


64 


4005,303 


309 


0.465 


0.193 


79016 


27443 


74 


4005.916 


998 


0.278 


0.084 


129218 


27205 


71 


4005,997 


18006 


0.248 


0.059 


910298 


215630 


21 


4005.1012 


2159 


0.230 


0.060 


564168 


147537 


27 


4005,1169 


1879 


0.521 


0.253 


555870 


230311 


53 


4007,2.4 


3 b i 


■'.300 


0.079 


64308 


1:410 






7 93 


-.324 


C.’Oi 


122553 


19596 




- 7 . 51 j 


7 " ' 
_ -• 


•). 333 


0. :50 


103O! 7 


33522 
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Tab Is C-4 Event Strength Suaaarv 



Nase 


r 


F'o 


si qua 


Fo 


sig&a 


(Tape 4, event #) 


U) 


(Pai 


•Fa) 


(N) 


IN) 


4009, 100 


3955 


0.217 


0.043 


237677 


47362 


4009,172 


370 


0.449 


0.240 


34708 


34086 


4009,178 


987 


0.256 


0.108 


150340 


52260 


4009,193 


2755 


0.234 


0.055 


261773 


57028 


4009,194 


1160 


0.321 


0.060 


191163 


28220 


4009.273 


12253 


0.267 


0.045 


339643 


56951 


4009,292 


16051 


0.243 


0.036 


315602 


46536 


4004,33) 


11046 


0.260 


0.039 


335321 


4970? 


4009.334 


979 


0.311 


0. 108 


159935 


50432 


4009,335 


3389 


0.256 


0.049 


322932 


62179 


4009,397 


3174 


0.223 


0.088 


250131 


98245 


4009,410 


492 


0.479 


0.179 


131094 


33345 


4009,447 


7887 


0.293 


0.048 


461103 


7573? 


4009,476 


333 


0.400 


0.290 


139730 


61426 


4009,493 


2 jpp 


0.173 


0.028 


199141 


31577 


4009,579 


3651 


0.343 


0.075 


373537 


81458 


4009,603 


1971 


0.21B 


0.074 


230200 


84135 


4009,620 


■3376 


0.237 


0.047 


315200 


61634 


: 009.o21 


a?3 


0.236 


0.065 


98243 


18204 


4009, cc7 


2367 


0.320 


0.195 


414294 


234724 


4009,725 


2032 


0.375 


0.132 


388341 


138491 


4009,736 


305 


0.981 


0.466 


22413? 


32413 


4009.S02 


3597 


0.746 


0.618 


1050885 


355095 


4009.317 


7934 


0.313 


0.089 


459710 


128412 


4009.823 


966 


0.179 


0.035 


105661 


22672 


4009,943 


3419 


0.277 


0.039 


323295 


46413 


4009.949 


9004 


0.244 


0.048 


294434 


57917 


4 009 ,56' 


7565 


0.339 


0.060 


*081*1 


7255a 


40'j9,1041 


1997 


0.224 


0.089 


221127 


70482 


4009,1160 


573 


0.425 


0.132 


137233 


31215 



69 

71 

53 

64 

65 

61 

60 

59 

66 

65 

49 

73 

5? 

69 

5 Q 

56 

58 

54 

56 

64 

53 

53 

58 

64 

64 

a 4 

64 

63 



-161- 



T=ble C-4 Event Strength Susmary 



NiSS 


r 


Po 


siqiaa 


Fo 


51 glad 


(Taos 4, event it) 


(a) 


(Pa) 


(Pa) 


(N) 


(N) 


4011.168 


6071 


0.198 


0.040 


234221 


46798 


4011,169 


9193 


0.222 


0.052 


276205 


64448 


4011,215 


403 


0.276 


0.079 


63100 


10454 


4011,241 


2997 


0.302 


0.123 


323546 


133350 


4011.250 


984 


0.353 


0. 162 


193201 


65473 


4011.259 


7903 


0.395 


0.055 


491277 


68119 


4011.278 


1159 


0.273 


0. 102 


153349 


53900 


-Oil .253 




0. 235 


0.042 


2393i0 


4I“0<> 


4011,307 


7800 


0.244 


o. too 


307771 


126135 


4011,344 


2132 


0.139 


0.041 


218914 


47390 


4011,352 


11454 


0.170 


0.033 


202390 


33875 


*01 1,335 


*086 


0.162 


0.036 


207993 


45867 


4011,532 


612 


0.316 


0. 120 


95187 


24435 


4011.541 


650 


0.183 


0.033 


53821 


14260 


4011.595 


367 


0.708 


0.209 


15750B 


32S76 


4011,615 


7633 


0.253 


0.056 


291097 


64418 


4011.516 


2958 


0.314 


0.416 


486903 


607412 


■tOi 1 .6a j 


3375 


0.217 


0.037 


227709 


33947 


*011.654 




0.377 


0.083 


3179" 


1412! 


4011,693 


t>790 


0.266 


0.051 


365550 


70710 


4011,721 


3 


0.245 


0.067 


52222 


17894 


4011,811 


974 


0.244 


0.056 


120519 


224-5 


4011.363 


1015 


0.166 


0.033 


187514 


*6708 


4011.355 


1043 


0.362 


0.067 


2*4063 


*6466 


*• 01. : 85 


10375 


0.238 


0.042 


31561c 


5 a Go 1 


4011,992 


2294 


0.292 


0.116 


320173 


117435 


4011,1046 


2000 


0.215 


0.055 


229653 


33118 


-Oil. 1 123 


730 


0 


0.082 


118230 


20570 


4011,1137 


3004 


0.267 


0.090 


299387 


985C-4 



65 

62 

67 

69 

5? 

62 

6& 

74 

61 

56 

65 

59 

66 

70 

64 

67 

46 

7 1 

60 

5a 

7 " 

f i. 

45 

31 

54 

5E 

ci 

38 

70 

65 
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Taoie C-4 Event Strength Summary 



Nase 


r 


Po 


sigma 


Fo 


si Qffia 


(Tape 1, svsnt i) 


(a) 


(Pai 


iPa) 


(Ni 


(N) 


4013,120 


5039 


0.232 


0.046 


249882 


49666 


4013,177 


1932 


0.224 


0.047 


246066 


59194 


4013,277 


537 


0.291 


0.065 


74728 


14232 


4013,317 


1326 


0.203 


0.037 


162763 


32443 


4013.354 


1712 


0.229 


0.045 


201933 


38367 


4013,373 


716 


0.354 


0.145 


136530 


31449 


4013,381 


3079 


0.313 


0.065 


325157 


66566 


41 13.403 


s03v 


0.30c 


0.068 


36636? 


31069 


4013,464 


2413 


0.240 


0.058 


306733 


70881 


4013.473 


664 


0.247 


0.072 


70133 


14240 


4013,518 


522 


0.362 


0.115 


107996 


■ 28307 


4013.661 


7801 


0.240 


0.078 


260305 


84007 


4013,664 


1567 


0.246 


0.075 


193674 


52565 


4013,694 


9314 


0.241 


0.058 


262115 


62356 


4013,723 


3458 


0.230 


0.053 


268779 


62309 


4013,776 


9750 


0.20? 


0.040 


403424 


77526 


4013,797 


7024 


0.203 


0.050 


250307 


59593 


4013,801 


422 


0.251 


0.077 


79381 


15347 


4013,866 


309 


1.157 


0.556 


219930 


82699 


4013,375 


324 


0.519 


0.247 


122952 


56241 


4013,950 


2536 


0.251 


0.045 


350436 


63028 


4013.979 


7943 


0.214 


0.042 


360405 


70664 


4013,1014 


2894 


0.180 


0.043 


284500 


69680 


4013.1018 


590 


0.245 


0.076 


87526 


16507 


4013.. 107 


335 


0.485 


0.134 


206443 


51643 


4013,1124 


12137 


0.290 


0.071 


316572 


77457 


4013.1166 


7055 


0.222 


0.063 


320351 


90645 


4013,1137 


982 


0.355 


0.103 


139420 


31011 


4013,1203 


8126 


0.194 


0.043 


234018 


5239<= 


1015.232 


434 


0.663 


0.223 


161500 


48620 


4019.178 


19953 


0.312 


0.054 


382727 


66870 


4021,81 


17603 


0.238 


0.049 


349159 


71641 


4027,28 


1068 


0.283 


0. 195 


145002 


51904 


4029 ,6t6 


10875 


0.172 


0.027 


24979? 


39740 


4C2 : ,39' 


*64 


0.162 


0.052 


3340',' 


7374 


4029.10:4 


9695 


0.21 3 


9.046 


■j-j79bQ 


c6->!d 


4031,823 


405 


0.246 


0.102 


66153 


36474 



7! 

58 

64 

58 

64 

61 

71 

J4 

53 

75 

58 

71 

64 

71 

64 

40 

64 

48 

67 

j-j 

51 

46 

46 

71 

. a 

*_■- 

71 

53 

8 0 

fc 4 

53 

64 

53 

64 

£*+ 

■Jv 

64 
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Table C-4 Event Strenatn Suinfflsrv 



Naae 


r 


Po 


Hape 1, event #) 


(in) 


(Pa) 


4033 . 490 


4896 


0.232 


40 W • 1 0 -4 d 


10359 


0.190 


4033,1053 


309 


0. 365 


4040,143 


7355 


0.392 


4040,385 


4268 


0.463 


4040.393 


12116 


0.438 


4 1*4 -7 , / 1 n 


15100 


0.-47 


4 a 40 . ^30 


18291 


0.491 


4040.371 


9545 


0.591 


4040,1157 


10323 


0.457 


4047,984 


499 


0.759 


4049,1064 


8435 


1.082 


4049,1066 


8221 


1.281 


4061,379 


4157 


0.374 


4U63,3o4 


2412 


0 . 3 18 




12033 


0.562 


4063.936 


1452 


0.290 


ZOO 1 . 5*3 


6170 


0.436 


2009 , 455 


-’’74 


0.470 


2009 . 455 


2464 


0.405 


2009.923 


sSl 


0.2“'? 


3001.11 


3922 


0.274 


3001,301 


3081 


0.214 


3047,169 


19233 


0.237 


3047,751 


1040 


0.280 






N = 






lean = 



5td jev = 
Max = 



51 0sTl3 


Fq 


sigma 


t 


(Pa) 


(N) 


!N> 


(Hz') 


0.048 


386811 


80418 


46 


0.027 


674418 


94838 


22 


0.120 


68367 


19654 


69 


0.055 


1119069 


156837 


27 


0.092 


690510 


121747 


51 


0.119 


848671 


229537 


40 


0.118 


7jo9’9'7 ! 


19958! 


4 * 


0.105 


962200 


204663 


40 


0.093 


1061264 


167076 


43 


0.108 


736054 


135924 


45 


0.252 


214017 


54162 


60 


0.135 


3389218 


424279 


2 5 


0.242 


4938915 


931151 


20 


0.353 


1152659 


107727 


■1C 

LJ 


0.056 


545947 


9&jiv 


*1» ; 


0.068 


2040664 


245937 


21 


0.027 


505945 


8 1 o58 


27 


0.065 


699233 


94401 




0.057 


1294184 


lu7S6l 


wU 


0.06-5 


973961 


156751 


79 


0.041 


193006 


6i3^4 


• U 


0.045 


772736 


123159 


27 


0.039 


395029 


7194a 


40 


0.037 


640540 


100804 


29 


0.112 


150767 


56491 


64 




151 








431151 








ujj-f It 








4938915 








334-00 
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Table D- 



Table D— 



APPENDI X D 



: Angles, Ranges and Times for Refractive 

Propagat i on Paths 



i: Spreading Loss Function, G(r) 
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Table D-l Angles, Ranges and Times for Refractive Propagation Paths 

Theta zero Theta 1 Theta 2 Theta 3 Theta 4 z 

Theta zero is the surface launch angle. 

Theta 1 is the angle at 80 meters. 

Theta 2 is the angle at the hydrophone (93m). 

Theta 3 is the angle at 254 meters. 

Theta 4 is the angle at 363 meters, 

z is the maximum depth of the ray. 

(See Figure 3-8) 



0.0s57 


0.039 


0.001 


93 


0.0658 


0.039 


0.004 


93 


0.0660 


0.040 


0.006 


93 


0.0670 


0.041 


0.013 


94 


0.0680 


0.043 


0.018 


96 


0.0690 


0.045 


0.021 


97 


0.0700 


0.046 


0.024 


98 


0.0710 


0.04B 


0.027 - 


99 


0.0720 


0.049 


0.030 


100 


0.0730 


0.051 


0.032 


102 


0.0740 


0.052 


0.034 


103 


0.0750 


0.053 


0.036 


104 


0.0760 


0.055 


0.038 


105 


0.0770 


0.056 


0.040 


107 


0.0780 


0.058 


0.042 


108 


0.0790 


0.059 


0.044 


109 


0.0800 


0.060 


0.046 


ill 


0.0810 


0.062 


0.047 


112 


0.0820 


0.063 


0.049 


113 


0.0330 


0.064 


0.051 


115 


0.0840 


0.065 


0.052 


116 


0.0850 


0.067 


0.054 


118 


0.0860 


0.068 


0.056 


119 


0.0870 


0.069 


0.057 


121 


0.0880 


0.070 


0.059 


122 


0.0890 


0.072 


0.060 


124 


0.0900 


0.073 


0.062 


125 


0.0910 


0.074 


0.063 


127 


0.0920 


0.075 


0.064 


128 


0.0930 


0.077 


0.066 


130 


0.0940 


0.078 


0.067 


131 


0.0950 


0.079 


0.06? 


133 


0.0960 


0.080 


0.070 


135 


0.0970 


0.081 


0.071- 


136 


0.0980 


0.083 


0.073 


138 


0.0990 


0.084 


0.074 


140 


0.1000 


0.085 


0.075 


141 


0. 1010 


0.086 


0.077 


143 



) 
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Angies tor Retractive Pr ocagat ioo Paths 
beta zero Theta i Theta 2 Theta 3 Theta 4 



0.1020 


0.087 


0.078 


145 


0. 1030 


0.089 


0.079 


146 


0.1040 


0.090 


0.081 


148 


0.1050 


0.091 


0.082 


150 


0.1060 


0.092 


0.083 


152 


0.1070 


0.093 


0.085 


154 


0.1030 


0.094 


0.086 


155 


0.1090 


0.095 


0.037 


157 


o. : :ov 


0.097 


3 . 088 


159 


0.1110 


0.093 


3.090 


161 


0.1120 


0.099 


0.091 


163 


0.1130 


0.100 


0.092 


165 


0.1140 


0.101 


0.093 


167 


0.1150 


0.102 


0.094 


169 


0.1160 


0.103 


0.096 


171 


0.1170 


0.104 


0.097 


173 


0.1180 


0.106 


0.093 


175 


0.1190 


0.107 


0.099 


177 


0.1200 


0.103 


0.100 


179 


). 1 Z !0 


0. 109 


0.102 


IS ! 


0.1220 


0.110 


0.103 


133 


0.1230 


0.111 


0.104 


185 


0.1240 


0.112 


0. 105 


137 


0. 1250 


0.113 


0.106 


189 


0.1260 


3.114 


0.103 


191 


0. 1170 


0 . 11 s 


0.109 


i 94 


0 . 1 2 E 0 


3.117 


0.1 10 


196 


,•.1290 


0.1 IS 


0.511 


198 


0. 1300 


■ 3.114 


0.112 


200 


0.1310 


0.120 


0.113 


202 


0.1320 


0.121 


0.113 


205 


0.1330 


0.122 


0.116 


207 


0.1340 


0.123 


0.117 


209 


0.1350 


0.124 


0.118 


212 


0.1360 


0.125 


0.119 


214 


0.1370 


0.127 


0.120 


216 


0. 1330 


0.123 


0.121 


219 


0.1390 


0. 129 


0.123 


221 


0.1400 


0.130 


0.124 


■i -37 

lL-j 


0 . 1410 


3.131 


0. 125 


21 6 


. 142 - 


0. 132 


o. : 2s 


T 23 


. ,-30 


0. 133 


0. 127 


231 


0.1440 


0.124 


G . 128 


233 


0.1450 


0.135 


0.129 


236 


O . li-iO 


0.136 


C . 130 


238 


' , J 70 


■ 3 . 


132 


241 
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Angles for Refractive Prooagation Paths 



Theca :ero 


Theta 1 


Theta 2 


Theta 3 


Theta 4 


0.1480 


0.138 


0.133 






0.1490 


0.13? 


0.134 






0.1500 


0.140 


0.135 






0.1510 


0.142 


0.136 






0.1520 


0.143 


0.137 






0.1521 


0. 143 


0.137 






0.1522 


0.143 


0.137 


0.002 




0.1523 


0.143 


0.138 


0.006 




0. .524 


1.143 


0. 138 


0.003 




0.1525 


0. 1^3 


0.138 


0.010 




0. 1526 


0.143 


0.133 


0.011 




0.1527 


0.143 


0.138 


0.013 




0.1523 


0.143 


0.133 


0.014 




0.1529 


0.144 


0.133 


0.015 




0. 1530 


0.144 


0.133 


0.016 




0.1540 


0.145 


0.13? 


0.024 




0.1550 


0.146 


0.140 


0.029 




0,1560 


0.147 


0.142 


0.034 




0. 1570 


0.143 


0. 143 


0.039 




'.1530 


0.14' 


0.144 


0.043 




0. 1590 


0.150 


0.145 


0.046 




0.1600 


0.151 


0.146 


0.050 




0.1610 


0.152 


0.147 


0.053 




0.1620 


0.153 


0.148 


0.056 




M630 


0.154 


0.14? 


0.05? 




). 1640 


0.155 


' . 150 


0.061 




0. la 50 


0.156 


0.151 


0.064 




0. IsaO 


0.157 


0.153 


0.067 




0. 1670 


0.139 


0.154 


0.069 




0.1630 


0.160 


0.155 


0.071 




0.1681 


0. leO 


0.155 


0.072 




3.1602 


0.160 


0.155 


0.072 




0.1683 


0.160 


0.155 


0.072 




0. 1684 


0. 160 


0.155 


0.072 




0.1685 


0.160 


0.155 


0.073 




0. 1686 


0.160 


0.155 


0.073 




0.1687 


0.160 


0.155 


0.073 




’■.1688 


0.160 


0.156 


0.073 


0.005 


0.1639 


0.161 


0.15s 


0.074 


0.003 


0. 1690 


0 . 1 s i 


0.156 


0.074 


0.010 


;.ii?i 


0. Sal 


0.156 


0.074 


0.011 


>. 1 58 1 


■; . 1 s 1 


v . 1 3 c 


0.074 


0.0*3 


0 .I 6 ?: 


0.161 


0.156 


0.074 


0.014 


0.1694 


0.161 


0.156 


0.075 


0.015 


0.1693 


0.161 


). 156 


0.075 


0.016 


’ It's 


0.161 


0.156 


9.075 


o.o :* 7 



243 

246 

248 

251 

254 

254 

25* 

255 

255 

256 

257 

257 

ICQ 

i-JU 

258 

259 

265 

272 

278 

284 

291 

297 

304 

310 

317 

323 

j u 

337 

343 

350 

357 

358 

358 

359 

360 

360 

361 

362 

364 

368 

371 
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Theta zero 
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0.133 
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0.2000 


0.193 


0.109 






0.2100 


0.203 


0.200 






0.2200 


0.214 


0.210 






0.2300 


0.224 


0.221 






0.2*00 


0.234 


0.231 






0. 2500 


0.244 


0.241 






0.2600 


>.255 


0.252 






0.2700 
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. . .75 


/. _ "7" 
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0.3000 
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ft. 3 18 
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395 

393 

401 

405 

408 

411 

415 

418 

422 

425 

428 

432 

435 

439 
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449 

452 

156 

159 

4*2 

46a 

489 

473 

507 

542 

576 

611 

647 

682 

718 

753 
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finales for Refractive Propagation Paths 


Theta :era 


Theta 1 Theta 2 


9.3600 

0.3800 

0.4000 

0.4500 

0.5000 

9.5500 

0.6000 

9.6500 

■,.'000 

0.3009 

0.9000 

1.0000 

1.1000 

1.2000 

1.3000 

1.4000 

1.5000 

1.5708 


0.356 0.354 
0.377 0.375 
0.397 0.395 
0.447 0.446 
0.497 0.496 
0.548 0.546 
0.598 0.597 
0.648 0.647 
'.633 O.o97 
0.799 0.798 
0.899 0.398 
0,999 0.999 
1.099 1.099 
1.199 1.199 
1.300 1.299 
1.400 1.400 
1.500 1.500 
1.571 1.571 
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Ranges for Retractive Propagation Paths 

Theta zero rl r2 r3 r4 r5 31 P2 

ri is the horizontal projection of the path frasi the source (Ohi) to 8Cm. 

r2 is the horizontal projection of the path fr om 30m to the hydrophone depth. 

r3 is the horizontal projection of the path from the hydrophone cepth to 254s. or the vertex 

r4 is the horizontal projection of the path front 254m to 362m, or the vertex, if sooner. 

r5 is the horizontal projection of the path from 362m to the vertex. 

Rl is the horizontal distance to the hydrophone on the downward path of the rav. 

R2 is the horizontal distance to the hydrophone on the upward swing of the r sy. 

*See Figure 3-8) 
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Retrsctive Propagation Paths 
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Ti*es for Refractive Propagation Paths 

tl l 2 t3 t4 t5 Ti 12 

tl, t2. t3, t4, and t5 correspond to the times required to traverse 
the oaths associated with rl, r2, r3, r4 and r5. 

Tl and T2 correspond to the times required to traverse the distances Rl and R2. 
1 See Figure 3-8) 
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Tiaes for Refractive Propagation Paths 
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Tiiriss tar Refractive Propagation Paths 
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Tiaes for Refractive Pronaqation Paths 
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Table D-2 Spreading Loss Function, 8(r) 

Theta zero Theta 1 z R1 R2 G(R1) G(R2) 10 log 10 log 

G(R1) G(R2) 

Theta zero is the surface launch angle. 

Theta 1 is the angle of the refractive path at the hydrophone depth (93m) 
z is the aaxiiua depth of the refractive path. 

R1 is the horizontal range to the hydrophone intersected on the downward swing 
R2 is the horizontal range to the hydrophone intersected on the upward swing 
S(R) is the spreading loss function. 
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Ssrsasinq losa Function, S(r) 
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